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COMBINED EXPANSIONS OF PRODUCTS OF SYMMETRIC POWER 
SUMS AND OF SUMS OF SYMMETRIC POWER PRODUCTS 
WITH APPLICATION TO SAMPLING! 


By Paut 8S. DwreEerR 


PREFACE 


This article is divided into two parts. Part I has for its title “Combined 
Expansions of Products of Symmetric Power Sums and of Sums of Symmetric 
Power Products” and develops the general mathematical theory which is ap- 
plied in Part II to “The Fundamentals of Sampling.’”’ Part II will appear in 
a latter issue of this journal. 

Each part is treated as an organic unit and has its own introduction and 
bibliography. Each article is assigned a given number and each book is given a 
letter so that references can be indicated concisely in the body of the dissertation. 

Each part is divided into chapters and sections. Braces are used to indicate 
the important formulas. 


PART I. COMBINED EXPANSIONS OF PRODUCTS OF SYMMETRIC 
POWER SUMS AND OF SUMS OF SYMMETRIC POWER PRODUCTS 


Introduction 


The mathematical material which is presented here has proved useful in 
generalizing that portion of the fundamental theory of sampling in which 
relations are established between the moments of the sample and the moments 
of the parent population. It is the purpose to establish the theorems in algebraic 
form since they constitute an extension of partition and symmetric function 
theory and may be of value to someone not necessarily interested in sampling. 

A great deal of work has been done in symmetric function theory but not 
much of this is of present value to the statistician. His problem deals with 
the “power sum’”’ while the classical theory, for the most part, deals with the 
interrelations of elementary symmetric functions and monomial symmetric 
functions. Only one phase of the reasoning developed in this investigation 
seems to have received extensive consideration previously and that is the subject 
covered in Chapter III. 

Previous authors have noted that much of symmetric function t!.eory re- 
duces, with a proper choice of notation, to partition theory. It is the plan of 
this treatise to present in Chapter I an outline of new partition theory which 


1A dissertation submitted in partial fulfillment of the requirements for the degree of 
Doctor of Philosophy in the University of Michigan. 
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2 PAUL 8S. DWYER 





shows how the parts of one partition are combined to form the parts of another 
partition, and which serves as a means of expressing the main result of Chapters 
II, III, IV, V. 

Chapter II shows how the formulas of Chapter I are applicable to the problem 
of finding products of power sums. The multiplication theorem for power 
sums, a generalization of the multinomial theorem, is stated in terms of power 
product sums and appropriate special cases are indicated. 

Chapter III deals with the expansion of power product sums in terms of 
power sums and shows how the formulas of Chapter I may be used. 

Chapter IV is the key chapter of the paper. The problem is to expand 
products of power sums in terms of power product sums, to multiply each 
power product sum by a quantity which is a uniquely defined function of the 
quantities composing the power product sum, and then to expand back in 
terms of all possible power sums. It is shown that the results can be written 
in a compact form which also utilizes the results of Chapter I. This result, 
as is shown in Part II, is directly applicable to the sampling problem of finding 
the moments of the sample moments in terms of the moments of the universe. 
Extensicn is made to multivariate distributions in Chapter V. 








Chapter I. The Combination of the Parts of Partitions 


It is the purpose of this chapter to provide a precise notation which shows 
how the parts of one partition of r may be combined to form the parts of another 
partition of r. For example, 2111, a four part partition of 5, can be made 
into 32, a two part partition of 5, by combining the three unit parts into a new 
part or by combining the 2 with one of the unit parts to form the 3 and the 
other two unit parts to form the 2. This last formation can be made in three 
different ways since anyone of the unit parts might be combined with the 2. 
The combination of the parts of the partition 2111 to form the parts of the 
partition 32 is to be indicated symbolically by Ps, + 3P22 where the subscripts 
indicate the number of parts collected and the coefficients indicate the number 
of ways in which an equivalent collection can be made. 






















1. Definitions and Notation. a. Partition [G; 105] [K; J; 1] [16; 105]. We 
consider the integer r to be composed of r unit indistinguishable parcels and 
define the partitions of r to be all those different groupings into new parcels, 
each new parcel containing one or more unit parcels, such that each resultant 
grouping of parcels contains exactly all the original r unit parcels. For example 
the partitions of 4 are 


4:31 ;22;211;1111 


b. Parts of Partitions. 
the parts of the partition. 


The numbers of the grouped unit parcels indicate 
Thus the partitions of 4 above 


4 ;31 ; 22 ; 211 ; 1111 have respectively 
3: 332s 8 3 





4 parts. 
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The pattern 22 may also appear as 2°. In general a p part partition of r is 
to be designated by 


Pi P2 Ps +++ Pp Where the p’s may or may not be equal and where p; + pz + ps 
4+.» +p, = ror by 


Di & Po *& Ps *& sR ++ RD 
- ps’ where< pir: + pom, + --- + Pome =7 
m++m+---+m=p 


c. Order of Partitions. When the parts of a partition are arranged in de- 
scending order we say that the partition is ordered. Thus 


¥1 


Pi 


Pipe --+ pp, is ordered if pi > po > ps > --- > Pp 


and pi: --- pr* is ordered if p: > pe > ps > --+ > Ds. 


For example, 21’ is an ordered partition while 312 is not. Unless otherwise 
- specified it is hereafter assumed that all partitions are ordered. 

It is sometimes convenient to refer to the order of the partition which is 
the size of the largest part, p, , when the partition is ordered. Thus the two 
part partition, 31, is of the order 3, while the four part partition, 1111, is of 
order 1. The set of the numbers pi' --- ps‘ is to be known as the complete 
order. 

These definitions of order and part are consistent with the usual definitions. 
[16; 105-106] [K; I; 1] [G; 100]. The concept of complete order, as far as I 
know, is not found in the literature. 

d. Weight of Partitions. Isobaric Partitions. The weight of any partition 
is defined to be the sum of all the parts of the partition. Thus the weight of 
Pi':+: pe’ iSPim1 + Pom2 +---+p.7, =7. Partitions having the same weight 
are called isobaric. Thus 4 and 211 are isobaric partitions. 

e. Algebraic Partitions. If the r original units are composed of a; , a, 43, 

- ad, nonseparable primary units, then the result of combining these in any 
possible way is to be called an algebraic partition since the r original units 
are now replaced by the r algebraic quantities a, , a2, --- a,. Thus aq, a2, as 
may be combined to form 


Gy + G2 + G3; a + 2-3; A + G3-A2; dz + 3-1; A1-d2-Gs 


‘which are the algebraic partitions of a; + a2 + a3. 

The parts of the algebraic partitions are the resulting combinations while 
the order and complete order, which indicate the numbers of algebraic expres- 
sions combined, agree with the order and complete order of the partitions in 
which the a’s are unity. The weight, which is equal to the sum of the parts, 
is indicated by w = a, + a2 + ----+a,. Thusif a, = 5, a = 4, anda; = 3, 
w=12. It is to be noted that the algebraic partitions are formed by combining 
the parts 5, 4, 3 and not by combining all parts of 12. 
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Now a; a2 --- @, is itself a partition of weight w = a, + a+ --- + 4,. 
If groups of the a’s are alike it may be written 










ai a2 


ay'az? --- a,” where 








Q;Q; + Aza, + +--+ + Aha, = Ww 


a + a+ +++ + a=. 











Algebraic partitions having the same weight are called isobaric. 


r 


f. Partition Combination Notation. Let Lee ' a indicate the number of 
are 


8 


different ways the r units, ordinary or algebraic, can be collected to form the 
5 


partition. Thus (32) indicates the number of ways in which the five units 










can be collected to form a partition with three units in one part and two units 
in the other. Since the three units forming the first part can be selected in 
sC3; ways and since this selection automatically indicates the other two units 
forming the second part, it follows that 


; 4 
ao = 5C3 = sC2'= 10. It is to be noted that LC 


yo) = 38 =6 


for if the four unit parts are a; , a2, a3, a4, then the three 22 partitions are 





Qi -+ Ae: Az + M4; G1 + G3-A2 + G4; Gy + 4-2 + Az 
since 








G3 + G4: + G2; Az + G4; + G3; Az + A3°Gi + Ay 


are essentially the same groupings as the first three indicated. 





2. Formula for is , a In establishing this formula we first take the 
Pt cel 
case in which no p rt is repeated. i.e. 7, = me --- 73 = 1 and py + po + ps 






--- + p, = r. In this case the formula becomes 


( 1’ ) = > 
Pip2-+-+ Ds Pi! po! +++ pe! 


‘Lis results from the fact that the p; units can be grouped in ,C>, different ways. 
The pe units then in ,_,,C>, different ways, the p; units then in ;_p,-»,C>p, dif- 
ferent ways etc. So that 


i 
, = Cp, ‘r—piC pg *r—p—paC ps py r—21—Pa---Pe-1C pe 
Pipe +++ Ps 


r! 


= >, Compare [B; 49][19; 12] 
Pi: pe! eee Ds* 





—_— vl «&$=eeto—lele 
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If however pi = po = --- = ps, then the same partition has been used s! 
different times since p, po, ---, Ps may be interchanged in s! different ways, 


so that 
Ores 
pi) (pr!)*s! 


1’ - etter etait its 
( oe ot) (pi!) (pal)? +++ (ps!)"* mil ma! ++ a6! 
Compare [19; 12, 13] [I; II, 252] 


By similar reasoning 


a;' a3? ay," 
‘1 2 a h 
3. Values of ( is oe Se 
Pi P2 +++ Ds 
ay: az* --- a,” may be collected to form the p parts of pf’ --- pz* may be indi- 


cated by 


) The number of ways in which the r parts of 


fay ag ah 
QA; Ag” +++ A, a 2111 fii 
® ™2 7 . I S aa = < 
Ur a... a :) Thus ( 39 ) 4 and ( 29 ) 3 
Formulas useful is evaluating this expression can be worked out from the results 


of this paper. A table of values of this expression for w < 8 has been given 
by the author [19; 29-32]. 


4. Notation for Combining the Parts of a Partition. Table I. We wish 
to indicate not only the number of ways in which a given r part partition of 
weight w can be grouped to form a p part partition of weight w, but also the 
number of parts of the r part partition grouped to form each of the p parts 
of the p part partition. As indicated in the opening paragraph, P3,; + 3P2. 


= ene, serves this purpose for the case in which the parts 2111 are col- 


ay ao ah 

aj! ag? +++ ap : . 

lected to form 32. P| 3, 7, ) serves this purpose in the more general 
Pi P2 eee Ds 

case. Its expansion gives sums of P functions whose subscripts are the numbers 

of parts combined and whose coefficients are the number of ways of forming 


the partitions from the parts. For example 


111 
P(*') = P,; ‘G) = 3P2; etc. 


e 
ah 


a) aq 
The use of P( “1%? “% ) is so fundamental to the present approach that a 
pi' pa? +++ ps" ' 


table is provided showing the different values when w < 6. 

Table I gives values of the function when w = 1, 2, 3, 4, 5,6. The values 
aj‘ --- a,” are given in the left hand columns and the values of pf' --- py* in 
the top row. The partitions are ordered from the top and from the left. To 
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, 2111 . 

find a given value, say P( 29 ) we note that w = 5, look for 2111 on the left 
and 32 at the top. The result is P3; + 3P:2. In the table the order of the 
subscripts is important in indicating the number of parts collected to form the 
respective parts of the ordered pj' --- ps°. 

The values in the table previously mentioned [19; 29-32] may be obtained 
when w S 6 by placing every P in Table I equal to unity. 

5: Value of P(ay' az’ --- a;"). The parts of the partition aj! az? --- ay” 
may be collected to form a large number of partitions of the type pr’ --- py". 


Thus the parts of the partition 2111 may be collected to form 5, 41, 32, 311, 
221, 2111. We denote by P(2111) the values 


2111 2111 2111 2111 2111 
r 4 ae 7 
w(2) +p) + ets eA) wa 
21th . i , 
+ P( 51,1) 2111 = Ps5 + 3Pu4l + [Pu + 3Px] 32 + 3Pou311 


+ 3Pe1221 + Pin2111 
and in general 


my 1 ah 
2 @;* +++ : 
P(af? af? --- af") = = P( 7 ) ll cca eli cca (2) 
Pi eee Ds 
where the summation holds for every partition pr‘ --- pJ* which can be formed 


by combining parts of aj’ --- a;". The values of P(af! .-- ay") for w S 6 
are given in the rows of Table I. Thus the value of P(2111) above is found 
along the row 2111 where w = 5. 


6. Values of P(1") and P(a’). When a, = 1 and aq; = r we have 


r ¥ TT T 
P(1’) = EP Yorn) +++ Ds’ 


1 


. 


. ] . . Tr rT 
and since there are ( ol different ways of forming p;' --- p;* from 
Ween 
the r units and each way is indicated by Py11...,r5 we have 


1 
P(1’) = cL 630g 3 
( ) mF. :) Pi py? Pl Pp { } 


When r = 2, 3, 4, etc., we get 

P(1’) = P22 + Py 1’ 

P(1*®) = P33 + 3P2 21 + Py 1° 

P(1*) = Py 4 + 4Py 31 + BP op 22 + 6Pey 211 + Pun 1° 
ete. 


as indicated in Table I. 





COMBINED EXPANSIONS 
Similarly when a; = a and a = r {2} becomes 


P(a’) = (0 ‘ a Pyjs...nte Capi)" (ape) --- Cap)" {3"} 


¥ ' ; . . 
since there are pe oo different ways of forming the partition (ap;)"' (ape)*? 
wee 


..+ (ap.)"* from the r equal a’s and each way is indicated by Ppr...pt+. 
For example 


P(a) = P(2)a = Pya 
aa 


2 2 
P(a’) = P(;,) + P(: ) = P,2a+ Pua 


3 
rad =#(S) +r (at) +4(¢5,) Pde ramen + Pa 


a-a-a 


P(a‘) = P,4a + 4P3; 3a-a > 3P 2 2a-2a oe 6P 2; 2a-a’ + Puna’. 


7. Values of P(a:a. --- a,). From the definition 


P(a;) = P() aq, = Py a, 
ay 


a, a2 a, de 
P(a;a2) = P(, + a + p(aes) 


= Pra, + a2 + Puasa, 


P(a,a2a3) = P( en Via + a&+ as) + p( sae \aF ara) 


a; + a3 + a3 Q1 + Ae- a3 


mas \ ( Q1a2ds 
+ P(. + =) , + a3-d2) + P — 


a1 42 43 


+ P(a1080) yas) = P;(a + de + a3) 


+ P| (ay + ae: as) + (a; + a3: as) + (ag+a3-a;)} + Piu(arazas) 


etc. 


Now if complete order of the general partition indicates the number of a’s 
collected to form the partition, the subscripts of the P’s are the respective 
complete orders. If we indicate the sum of partitions having the safne com- 
plete order by tie term “partition type” and indicate the partition type com- 
posed of all terms having the same complete order 


v1 


T a ad 
Pi . Ds . by i pil---pe! 
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P(a) = PiT; 

P(ai;a2) = P2T2 + PuTu 

P(a,az2a3) = P3T3 + PauTa + Pin Tin 

P(a,az2a304) = PgT,4 + PuTs + PoT2 + PouToun + PunTun 
etc. 


and in general 
Playa, «++ a) = Dy Pott...gte Tyit...p7e {4} 


This formula can be used in writing the formula of Table II or formulas of 
weight greater than 6. Thus 


P(543) is given by P373 + PaTa + Pin Tin 


where 
Ts = 12, T 2 = 9.3 os 8.4 t 7-5, and Tw = 504-6 


where the dots do not indicate multiplication, but merely the separation of 
the parts. 


In general 7',71...,z. is composed of * ' ol partitions since this is the 
2. ee 


number of ways in which the a’s can be combined to form partitions having 
the same complete order, pi' --- p;'. 
Formula {3’} is a special case of this formula. If the a’s are all equal, the 


~ 


—_ 1 
. o artitions are equal so that T7,2...,..=|[ ,; x, }(ap,)"! 
pit... pr p q pil---p5 rr: Pi) 


--+ (ap,)"*. Substitution in {4} gives {3’}. Similarly {4} gives {3} when 
all the a’s are unity. 


8. Generalization from Symmetry. The function P(aaz --- a,) is a sym- 
metric function of the parts a; , a2, --- a,, i.e., the interchange of any two of 
the parts does not change the value of the function. It is possible to use this 
fact as a basis of generalization and to derive {4} from {3} by its use. From 
{3} we have 


i 
P(1’) = . w, |) Pott...9%¢Di' +++ De" 3 
(1’) Gh. «al pilesple Di ++ D {3} 


where bas , a is the number of the equivalent partitions which can be 
—e 


formed from the r units. In case the r units are replaced by the r different 


+ 


a’s, there will result Le . “ different partitions having the same complete 


1 ove 3 
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order. These Ln ; oa different partitions defined by T'p1...,12: replace the 
oof 
we i pt equivalent partitions ot {2} and we have 
Pi eee 8 


P(aiaz --+ a,) = 2. PF oft..<ghe T gft...g8e {4} 


9. The Recursion Rule. It is possible to establish a recursion property by 


which the value of P(aia2 --- a,a,41) can be obtained from the value of P(a, a, 
..-a,). We note, from the results of Table I or by {4} that 
P(3) = P,3 


P(321) => P36 + P51 a P42 ~ P3233 oe Pyy,111 


P(32) is ubtained from P(3) by symbolic multiplication of its expansion, P:(3), 
by the expansion of P(2), P,(2) This symbolic multiplication ‘s s:co-.nvlished 
by adding the 2 to the 3 and also suffixing the 2 to the 3. If the 2 is added, the 
subscripts of the P’s are added while if suffixed, the subscripts of the P’s are 
euffixed. 

More generally if P(a;) = P,(a;) and P(az) = P,(a2), then the result P(a,a,) 
= P,(a,; + a2) + Pu(a:) (a2) is obtained by multiplying P,(a;) by P2(a2) [or 
P2(a2) by P:(a,)] symbolically if the subscripts are added when the a’s are 
added and suffixed when the a’s are suffixed. Similarly P(a,a2) = P2(a,; + ag) 
+ Py:(a;) (a2) when multiplied by P(a3) = P;(as3) gives 





P(a, 424s) = P;(a, + de + as) + Poa; + Ga° Gs) -+ Po(a, + G3: dz) 
+ Pr2(ai-a2 + as) + Pin(aia2as) 


when the rule of multiplication is the adding of a3 in turn to every part of every 
partition with the appropriate adding of subscripts and the suffixing of a; to 
every partition with the corresponding suffixing of subscripts. It is important 
to note that the P coefficient of a;-a_ + a3 is Py and not P2; although the term 
could be written Pea, + a3-d,. The applications do not demand the retention 
of a given order of subscripts though the continued application of the recursion 
rule does demand it. 

In general the value of P(a, --- a,a,41) can be obtained from the value of 
P(a,a2 --- a,) by the symbolic multiplication of the expansion P(a,a2 --- a,-) 
by P,(a,+1) since all possible algebraic partitions of a; + a2 + --+ + G + Gra 
are obtained from all possible algebraic partitions of a; + a2 + --- + a, by 
adding a,,; in turn to each part of each partition and by suffixing it to each 
partition. The corresponding P subscript, indicating the number of a’s col- 
lected, is increased by 1. 

The recursion rule is useful in checking the entries of Table I. As a matter 
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of fact Table I was computed with its use and the order of the subscripts is 
that which results from its use. The rule is also useful in finding values when 
w> 6. For example, since 


P(321) P36 + P51 -+ Py 42 + P33 + Pin321 

P(3221) = Ps8 + P3162 + P71 + P2253 + Pen512 + P3162 + P44 
+ Po 422 + P2253 + P3385 + P1332 + Pen521 + P1341 + P2323 
+ P3212 = P88 + P3171 +2 P62 + (Psu + 2P22)53 + Po44 
+ 2P21521 4+ Po1431 + Pon422 + 2Pon332 + Pi13221. 





A useful check is based on the fact that the sum of the P coefficients of 
P(a, --- a,) should equal the sum of the coefficients of P(1’). In the above 
illustration the sum of the coefficients is Py + 4P3, + 3Pe2 + 6Pea. + Pim as 
desired. 





10. Use of the P Function Formulas. The P function formulas, as defined, 
represent concisely the ways in which the parts of a given partition may be 
combined to get the parts of other partitions. They are also useful in writing 
expansions of certain partition functions whose expanded values are expressed 
in terms of other partition functions. They are used, in this paper, in expressing 
the multinomial theorem, the multiplication theorem for power sums, the 
expansions of power product sums in terms of power sums, expansions of mono- 
mial symmetric functions in terms of power sums, the double expai:sion theorem 
itself, the coefficients in the double expansion theorem as well as the sampling 
laws of Part II. They are also useful in representing the expansions of different 
moment functions and can be associated with important concepts of mathe- 
matics and statistics such as, for example, the differences of 0. Such applica- 
tions, however, are not pertinent to the line of reasoning which is developed in 
Chapters IT, III, IV, V. 


Chapter II 


It is the purpose of this chapter to obtain formulas for the expansion of 
power sums. 


11. Definitions. a. Power Sum. Let zx be a variable which is restricted to 
the N variates, 21, %2, 23, --- , Zw. Then the a-th power sum of the variable 
indicated by (a) is defined to be 





” 

(a) =a t+ait---+ay= Dat - {5} 

It is assumed for the purposes of this paper that a is a positive integer or 0. 

b. Power Product Sum. The expression )_ 27'2%? is to be called a power 
ee 
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product sum since it is composed of the sum of products of the powers of the 
variates. It is to be denoted by (a:-a2) or (a,a2). Thus >, 2327 = (3-2) or 
ii 


(32). The value (a-a) = (a’) = x 2{2z; is a special case of (a,a2) where 
ty 

a2 = a, = a. In general the power product sum is defined by the right hand 

member and indicated by the left hand member of 


(mm ---«@) = } ata -.- ae {6} 
£, ig igh- +s Mi, 


If 7; = %, the power product sum becomes 


(a; + de-a3-a4 +--+ a,) = » Uji Vig +++ Xe {7} 
ipmigheigheighe: + Mi, 

There are many different definitions since there are many different ways of 

indicating equality relations among the 7’s. Each results in a unique power 

product sum which is to be called, for brevity, a power product. If the a’s 

are all unity, there are many duplicates. Thus for the grouping p;' --- pry’, 


there are gt gle equal power products (pj! --- py‘). Inthe more general 
‘ee 


case we can let 7',71...,%. represent the pi? : al different power products 
PP ows 


8 
Te 


having the same complete order, pi‘ --- p; 
these forms having this complete order by 


We may represent any one of 


(919293 --+ Q%) 


where Q+e2+oet-::-+%=w 
. or by (qi'g2” --- gt‘)* 

where Mi%1 + Qetet+ +--+ +H = w 

and M+ eee tH = Op. 


c. Symmetric Functions. Both the power sum and the power product are 
symmetric functions of the variates since the interchange of any z; with any 2; 
does not change the value of the function. Also the powder product having 
p parts is composed of N® products of powers since the first group of equal 
z's may be selected in N ways, the next group in N — 1 ways etc. 


d. Monomial Symmetric Function. It is customary to use the monomial 
symmetric function which is defined as 


Qi 92 Ip 
snl i _ Vi, Vig ++ Ti, 
ti<te<ct3<"**<ty 


* “Tt was intended that the letter representing the exponents of the q’s should be the 
Greek ‘chi,’ and not the English ‘z.’’’ 
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and which we designate by M(q --- q,) or by M(qi' --- gi‘). This function 
is not useful for our purposes since the number of terms in its expansion varies 
with the number of repeated q’s. For example if N = 3 and q & q2; M (qq) 
= afta + atlas? + xia? + agtal? + xa? + zs' ae? = (gm) 

while if g; = g2 = 

) 


2 
2 
(q) = xiat + zizi + x23 = “or 


The monomial symmetric function keeps the number of product terms a 
minimum by eliminating all repeated terms while the power product sum 
keeps the number of product terms the same by the use of repeated terms, 
when some of the parts are alike. 


12. The Formula Connecting (qi'q2? --- qi!) and M(qi' --- qi‘). The 
power product is composed of N® products, each of which is repeated x!2»! 
N” 

Xi! Xo! +--+ 2,! 
different products which, when repeated 2!2! --- z,! times, gives the N” terms 
above. Hence 


-++2,! times. The monomial symmetric function is composed of the 


(qi! «++ gi) = aitae! --- a! M(qi' --- gi) {8} 
Zz zs 1 z1 ze | 
M(qi ee eee +++ gi’) {9} 
In the special case in which gq, = 1 and 1, = p 
e 
(1) = p!M() and M(1’) = © £10} 


The function, M(1°) is commonly called an elementary symmetric function. 
We refer to the corresponding (1°) as the unitary power product sum. 





13. Correspondence of Partitions and Power Products. To each power 
product (qi' --- qt‘) there corresponds an algebraic partition gj’ --- gi‘ having 
p parts and weight w = a, + a+ --- +4,. 

This follows at once from the definitions and notation. Thus if w = q 
+ a2 + a;, the power product 


Dd mena = QL anh att = (a + aa) 
iy=ig iz i ie 

is, by notation, associated with the partition a, + a2:a3;. Conversely each 
algebraic partition, when enclosed in parentheses, represents a power product 
sum. 

This proposition is useful in that it enables one to establish a relationship 
between the theory of power product sums and the theory of partitions. Earlier 
writers have used a similar correspondence in relating the theory of monomial 
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symmetric functions to that of partitions. See for instance [3; 106], [4; 5] 
(5; 1; 7]. 

Due to this correspondence we do not hesitate to apply such terms as part, 
order, complete order, similar, etc. to the power product as well as to the 
algebraic partition. Also the sum of all power products (gj' --- qi‘) having 
the same complete order is represented by T(pi' --- pz‘). This represents 


r 


the sum of * . wel similar power products. 
Pe os 


14. The Multiplication Theorem for Power Sums. The correspondence 
property enables us to derive a theorem, to be known as the multiplication 
theorem, which expresses products of power sums in terms of power products. 
The type of argument is introduced by establishing simple cases of the theorem 

= os a2 a 
(a;)(a2) (2 zi \(2 tj; ) ~_ du at aa? 
i1= 
to=l 


= Do aia + 2 aiay = (ai + a2) + (aie) 


t1=t2 iia 


(aaa) =(EaP(Ea")\(CaP) = YC attattals 


t1,429.%3 
= (a; + a2 + a3) + (ai + ae-a3) + (a; + a3: a2)’ + (a1: a2 + a3) + (a1-a2-a3) 


since the value p is broken into 


t1,t2,%3 
es a a % a 
t1=iomi3 i1=ig is i y=i3 ig i; Sig=is i, ig iz 
In general, when r S N 


(a)(m)--- (a) = 2) aia --- of 


a ..m 
and this can be broken into summations featuring different equality relations. 
These summations define all the different power product sums of weight w = a; 
+a2+---+a,. The different algebraic partitions of a; + a2 + a3+ --- +4, 
correspond to the different power product sums. It follows at once that the 


value of (a;)(a2) --- (a,) is obtained by writing each algebraic partition of 
a, + a + --- + ,, enclosing it in parentheses to represent a power product, 
and adding. More symbolically we have 
° (a;)(a2) --- (ar) = D> (gi! --- gi') {11} 
where qj' --- gt‘ represents any algebraic partition of a; + a + --- + a. 
and the summation holds for all such partitions or by 
(a:)(a) --- (ar) = 20 TQp --- pi") {12} 


~ 


where T(pi' --- pz‘) represents the| , : x, ) Similar power products and 
Pi © eid p 


the summation holds for each different complete order. 
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For exauple (a;)(a2)(a3) = T(3) + 7T(21) + T(111) 
and T(3) = (a; + az + as), T(21) = (a; + ag-a3) + (a; + a3-a2) + (a2 + a3-01), 
and 7T(111) = (a;-ae-as). 

The theorem has been established on the assumption that r < N. If such 
is not the case it is possible to satisfy the assumption by adding additional 
variates, tvi1, tvi2, --- 2, all 0, without changing the value of the power 


’ sums or of the product of the power sums since the added terms are always 0. 
Thus 


(xt + 22)(xt + a2)(ai + 22) = (at + x2 + xs)(2t + 22 + 25)(2i + 22 + 25) 








when 2; = 0 
Then 


(a)\b)c) = @+b+c)+ @+b-c) + @+c-b) + (+ c-a) + (a-b-c) 


which is 
b b 
Do asthe 4 De atte 4+ Datta + Dates + DS xtab at 
ij is iti ixp&k 


The term 2 x3} x, = 0 since every product composing it contains an z; = 0. 
ip Rk 
The other power product sums are to be applied to the original variates only 
since the terms involving 2; are 0 in every case. 
In general, if r > N, it is only necessary to write out the power product 
sums having WN or less parts since all those having more than N parts will be 0. 


15. The Multiplication Theorem Using the Results of Chapter I. Com- 
parison of {12} with. {4} shows that {12} can be obtained from {4} by placing 
P(ay --- a) = (a:)(a2) «++ (ay), Typs...te = TORT! +++ pl") and Py...gte = 1. 
Since this can be done for all values of a and r it follows at once that the entire 
theory of Chapter I is applicable to the present problem. For example Table I 
shows that 


P(821) = P36 a P51 ae P42 + P3233 - P3321 
and it follows that 
(3)(2)Q.) = (6) + (1) + (42) + (3) + (821) 


It should be noted that it is possible to use the table previously published 
{19; 29-32] since the entries in’ this table are the values obtained when 
P,11...p% = 1. The value (3)(2)(1) may also be checked from this table. 





16. The Multinomial Theorem. The multinomial theorem is a special case 
of the multiplication theorem for power sums in which the power sums are all 
equal. Ifa, =a@=--- =a, =1, 


® Be 1’ ® Ts 
T(pi' --- ps‘) = (peo ose )OOr “++ ps‘) 
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and {12} becomes 
OY =X (pn gee or BF (13} 


which is the multinomial theorem in terms of power product sums. Special 
cases are 


(1)’ = (2) + (41) 

(1)* = (8) + 3(21) + (11) 

(1)* = (4) + 4(31) + 3(22) + 6(211) + (1111) 
etc. 


The result of {13} may also be obtained immediately from {3} by placing 
P(1’) = ()’, Pon...pre = 1, and pr' --- prt = (pr* --+ pr). 

A more general form of the multinomial theorem is that in which a; = a 
=.--- =a, =a. In this case 


(xi + a2 + --- + 2h)" = (a)’ 
and {12} gives 


CY = Epes 2 pps )Mend" => Con") (14) 
where ((ap;)"! --- (ap.)"*) has parts api, --- ,ap.. Thus 
(a)* = (8a) + 3(2a-a) + (a’) 
so that 


(2)* = (6) + 3(4-2) + (2°). 


The result {14} may also be obtained immediately from {3’} by placing P(a’) 
= (a)', Pops... = 1, and (ap,)" --- (ap,)"* = ((ap))" --- (ap.)"). » When 
N = 2, {13} gives the binomial theorem 


ay = (0 Now 
special cases of which are 
(1)’ = (2) + G1) 
(1) = (3) + 3(21) 
(1)° = (4) + 481) + 3(22) 
(1)° = (5) + 5(41) + 16(82) 


etc. 
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. 
These can be readily translated to the usual form. Thus 
(a + b)* = a + Bf + 4(a°b + B’a®) + 3(a°b’ + Ba’). 
In a similar manner the trinomial theorem appears as 
ig 
ly = 
gy => Ls >, (Papers) 


A special case of the multinomial theorem {13} is also useful in writing N’ 
in terms of sums of N’. When the variates are all unity the power sums are 
all N, and the power product sums are the number of terms in the partition 


representing it. If a partition has p parts the number of terms in it is N™, 
We then have 


a, 
N’ = ( i ..) N° 15 
p> Pi -:: ps’ { 
Special cases are 
N aus N + N® 
N N + 3N + N® 
Nt=N+4N” + 3N® + 6N® + N® = N+ 7N® + 6N® + N™ 


etc. 





















17. The Use of Monomial Symmetric Functions. It is possible to express 


the results in terms of the monomial symmetric functions by means of {8}. 
Thus 


(2)(2)(2) = (6) + 3(42) + (222) 
= M(6) + 3M(42) + 6M(222). 


In general, Table I may be used to express products of power sums in terms 
of monomial symmetric functions. It is only necessary to place every Ppr1...p% 
= 1 and to multiply by the factorials indicating the repeated entries at the 
head of each column. The table [19; 29-32] may be used similarly. 

The multinomial theorem in terms of monomial symmetric functions becomes 


. Zi ® . 
CY = (seer gee) eateal oe «BP 
Pi eee Ds 
and by {1} ; 









_— a. oe a rare a, 
Y= 2 Gyagape = pale MOP + Ph) (3) 


as it is conventionally stated. 
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18. The Multiplication Theorem from the Multinomial Theorem. It is 
possible to use generalization from symmetry in deriving the multiplication 
theorem from the multinomial theorem though this can not well be done from 
its conventional statement (16). The monomial symmetric function: does not 
have the property that M(a-b) = M(a-a) when b = a while (a-b) does become 
(a-a) when b = a. The first step then is to reduce {16} to power product 
sums by means of {9}. We then have 


r! 
1)’ = compnenininiiaieienmniemnn 71 te 
( ) Zz (pi!)** (pe!) *2 age (p.')** a1! 39! eo r,! (pi Ps ) 
Next it is necessary to introduce the factor Le ’ el for there are many. 
oss wil 


equal terms for each value pr’ --- py* when the a’s are all unity. This is very 


® 


easy in this case since the value of the coefficient of (pf! - - - pi*) is DD ' wl 
Pe... ged, 


It follows at once that 


Qa) = 2 Le 7 ot)" +++ pr’). 


71 


Ri *** De 


Suppose that the r units are replaced by aia, --- a,. Then the ( ; “ll 
power products, (pi' --- pz*) will be replaced by the * ; al different 
f... ef 


power products composing T(pi' --- pz‘). It follows at once that 
(a:)(a) --- (a) = DY TQpi' --- pe). 


19. The Determination of the Coefficient of a Given Power Product in the 
Expansion of a Product of Power Sums. In some cases we wish to determine 
the coefficient of a given power product without computing the complete 
ra eee ah 

expansion. This is given by P a 2) where the P coefficients are unity. 
al 

Thus the coefficient of (32) in the expansion of (2)(1)(1)(1) is found from 


P cy = Pu + 3Pq and is 4. 


20. Relation to Previous Results. The multiplication theorem may be 
viewed as a generalization of the multinomial theorem. A more general proof, 
applicable to multivariate problems, could be presented with the use of more 
involved notation. It seems wise rather to present the simpler one variate 
case and to emphasize the principle of generalization from symmetry which 
will enable us to write the multivariate laws with relative ease. 

The general problem discussed here seems to have received a very small 
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amount of consideration as much of the extensive classical theory of sym- 
metric functions is limited to the interrelations of the elementary symmetric 
functions and the monomial symmetric functions. 

A monumental work on symmetric functions not subject to this limitation 
is the Combinatory Analysis of MacMahon [K]. MacMahon provided a 
technique for multiplying power sums in many variables as a special case of a 
more general theory. [K; II, 321]. 

Some of the work on alternants is closely related to the problem of products 
of power sums although the alternant, as usually defined, is limited to the case 
in which r = N [I; II, 446]. For an example the reader is referred to a devel- 
opment by Muir [L; 335-6]. 

Thiele (1889) gave tables’ of products of power sums in terms of monomial 
symmetric functions for partition products of weight £8 [H; 114-117]. J. R. 
Roe has later given one for w S$ 10[N; Plates 17, 18]. Statisticians have some- 
times stated the results in nontabular form. See for example, the multipli- 
cation formulas of Church [13; 81-83] [14; 370-1], whose results may not at 
first appear to agree with those above since Church has used a less compact 
notation and, of course, the monomial symmetric function. 

The chief contributions of the present attack are 

1. The use of the formulas and tables of Chapter I in writing expansions of 
products of power sums. 

2. The.use of power product sums in place of monomial symmetric functions 
which makes feasible. 

3. Generalization from symmetry. 




























Chapter III 


It is the purpose of this chapter to establish formulas giving the expansion 
of power products in terms of products of power sums. 





21. The Binet (Waring) Identities. It is customary to introduce this subject 
with formulas for M(a-b), M(a-b-c), etc. so we first derive the formulas for 
(a-b), (a-b-c), etc. We may use the results of Chapter II since the problem 
here is the inverse of the multiplication problem. By the multiplication 


theorem 
(a)(o) = (a + b) + (a-b) 
(a)(b)(c) = (2+ 6+) + (+ b-c) + (a + c-b) + (6 + c-a) + (a-b-c) 
(a + b)(c) = (@+b+c) + (a + b-c) 


2 These tables are not accessible to me, but Thiele refers to them in his ‘“Theory of 
Observations.”’ 
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so we get 
(a-b) = (a)(b) — (a + db) 
(a-b-c) = (a)(b)(c) — (@ + b)(c) — (a + c)(b) — ( + c)(@) 
+ 24a + b+ ¢) 

Similarly 
(a-b-c-d) = (a)(b)(c)(d) — (a + b)()@) — (a + c)()@) 

— (a + d)b)(c) — (6 + )@)@) — 6 + d@l©) 

— (€ + d)(a)(b) — (a + b)(¢ + a) — (@ +c) + dd) 

— (a + d)Q +c) + 24+ b + c)(d) + 24 + b + dc) 

+ 2(a + c + d)(b) + 206 + ¢ + d)(a) — 62 +b+4+c¢+4d) 

eS ah en eT ee ee {19} 
When a *& b *& c & d, {18}. {19}, {20} are also the formulas for M(ab), M(abc), 
M(abed). ‘These formulas are quite commonly attributed to Binet who gave 
them in 1812 in connection with certain proofs of determinant theory [1; 284] 
[I; I; 81]. Waring should be given credit (see Miscellanea Analytica 1762). 
Binet gave no proof. The reader is also referred to the earlier work of Paoli 
[A; section 28]. 

A much more adequate treatment was given by Hirsch in the early 19th 
century [B; 35-38]. He wrote, out the terms for M(a-b-c-d-e) and indicated 
a scheme for extending the results. More than this he proved that any “numeri- 
cal expression’’—his term for monomial symmetric function—can be reduced 
to numerical expression having one less part [B; 26]. The continued application 
of this theorem leads eventually to numerical expressions having only one part, 
i.e. to power sums. Hence ail numerical expressions can be reduced to power 
sums [B; 27, 32]. 

Recent authors give essentially the same proof. See for example Bocher 
[J; 241-242] who states the theorem, “Every symmetric polynomial is a linear 
combination with constant coefficients of a certain number of the 2’s.”’ See 
also O’Toole [16; 114] and Burnside and Panton [E; 167]. Thus modern authors 
provide a proof of the fact that M(a; --- a,) can be expanded in terms of power 
sums but most of them fail to provide a formula giving this precise expansion. 
Even MacMahon after writing the values of M(Au), M(dur), M(d’), M(a*) 
avoids the immediate generalization by stating [K; I; 7], “In actual practice 
there are easier ways of calculating the many part functions and the general 
formula is of little importance.” 

While MacMahon’s statement has a certain amount of truth in that any 
given monomial symmetric function may be computed from others having one 


less part by the recursion property described by Hirsch, yet there are many 
cases in which a definite formula, rather than a method, is desirable. A formula 
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particularly is demanded by the statistician who is working with a large number 
of monomial symmetric functions simultaneously. See for example the remarks 
and efforts of Carver [15; 103-104, 119-120], Church [14; 373, 377-378], and 
O’Toole [16; 115]. 

Some authors have provided solutions and it appears that statisticians are 
not entirely familiar with all the work whick has previously been done. It 
is the aim of the remainder of this chapter to suggest references which make 
previous work available to statisticians as well as to present a logical and quite 
complete development. The main results are not essentially new although 
their explicit statement in the language of power products is necessary for the 
development of the next chapter. The argument features the easy generaliza- 
tion from symmetry. The value of (1’) is expressed in such a form that the 
value (a; --- a,) may be obtained immediately from it. 






22. The Value of (1") from Waring’s Expansion for the Elementary Sym- 
metric Function. We first derive the formula (1’) from the conventional 
Waring’s expression for p,, in terms of the power sums. Burnside and Panton 
[E; II; 92] give this as 





rytrot-+-:t+rm 71 72 


ow (—1) 81" 83" +++ Sn" 


Pr 1) P21) > Fm f 1) 237 me om 


where pm = (—1)”M(1”) and where Sj! --- Sj is any 7) + r2 + --- + Tm 
part partition of m. When m =r and p;' --- p;‘ is any p part partition of m, 
{20} becomes 


(—1)’M(1’) =n by (—1)?(p)"! poacee (p.)*' {21} 


pil™! eee ps!™* ay! re! 7a T,! 


Dividing by (—1)’ and noting that (—1)”” = (—1)"” we have 


{22} 


May = DV" 


+ pol™* mi! m2! +--+ 1,! 


and hence that 





(1’) i r!M(1") = 2 (—1)"’r! (p:)" oe (p.)** {23} 


pil™! +++ py! **ay! 29! - ++ 2! 


A second proof of {23}, given in the next sections, does not assume the 
formula {20} and develops by easy stages. Although somewhat longer than 
the method above, it contacts much of the work that has been done in this 
field. It also provides two useful arithmetic checks dealing with the coeffi- 
cients which the more analytic method above does not provide. Those who are 
familar with {20} above and are interested in the immediate development of the 
argument with the use of {23} should turn to the equivalent {38} of section 28. 
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23. The Newtonian Formulas. The development begins with the well 
known formulas connecting the power sums and the elementary symmetric 
functions which appeared in Newton’s Arithmetica Universalis. These formu- 
las are given by Bocher (J; 244) as follows 


Si — PiSea +--+: + (—1)*" pet St + (—1)* kp; =0 kK=1,2,--- {24} 


where S; is the sum of the k-th powers and p is the 7-th elementary symmetric 
function. 

So many proofs of this theorem are accessible that a repetition here is hardly 
justifiable. A proof using calculus was given by Bocher (J; 243). Proofs 
using algebra only were given by Hirsch (%; 16) and Chrystal (F; I, 437). 
Muirhead (9; 66-70) gave three proofs of » ‘1ich the second is perhaps best 
adapted to the present development. 


24. The Determinant Equivalent of (1’). It is usual to solve the Newtonian 
equations for the power sums (J; 244) but our objective is the solution in terms 
of the power sums. The equations are 


Pi 

(1)p. — 2p2 

(2)p2 — (1)p2 + 3ps 

(3)p: — (2)p2 + C1)ps — 4p = (4) 


whence 


0 (1) 
«fj 0 (2) 
—(1) (3) 


(-1)""~r -1) @-)) 


(r — 2) -(@ — 3) G — 4) 
_ |@-1)- @ - 2) @ - 3) Sa 9) (r)| 
1 0 0 0 


(1) —2 0 
(2) —(1) 3 
(r — 2) — (r — 3) (r — 4) (—1)"(r — 1) 0 
(r — 1) — (r — 2) (r - 3) (—1)"*(1) (—1)""r 


Next, factor out all the negative signs in the even numbered columns in each 
determinant. The number of these columns of negative signs is the same as 
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the number in the denominator if r is odd. If r is even, there is one more in 
the denominator. - Hence the negative signs may be dropped in both deter- 
minants if (—1)"" is inserted in the numerator. Furthermore the value of 
the determinant in the denominator is r! Next, change the numerator by 
moving the r-th column to the first column position and inserting the com- 
pensating factor (—1)". If A, represents the resulting numerator determinant, 
the value of p, becomes 


», = (D=D" 





r! ~ 
and 
A, = r!p, = (1) 
We have then 
|) 1 Me eee eee 0 0 | 
@)  @) ee 0 0 
ae eee Tl 
a eines Q) r—1) 
| (r) a) a (2) (1) | 


The determinant has received the attention of earlier writers {19; 3}. Gen- 
eralizations of it will be mentioned at the close of the chapter. Its expansion 
in terms of power sums is known and may be written 





_ yy (-1)?r! (p)” --- a” 
tin ee eae 126} 
where Pit a P2 Te +.--- + DsT, = 7 


and ™m + m2 +--+: +, =p. 


See for example O’Toole (16; 113). 
It is at once evident that {26} is equivalent to {23}. Those who are familiar 
with the expansion of A, above may wish to turn immediately to {38} of sec- 
tion 28 since the intervening sections are devoted to a rather detailed and rigor- 
ous expansion of the determinant. This development follows, in a general way, 
that given by Mola (5; 190-195). 


25. The Expansion of (1°) = A,. The determinant, A,, is a special type 
of determinant which is known as a recurrent. There is a simple recursion 
property which is useful in its expansions in terms of products of power sums. 
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(1) 1 Oh hicehitaiaiiis 0 0 
(2) ~=) sdinnsdoeeieacs 0 0 
icine (3) (2) Ee | °ditslaaad aibabaeain 0 0 | 
 @-)@-2) GQ) or 
(r+ 1) (@) (r — 1) 2) @Q) 
If we expand A,,,; in terms of the (r + 1)st column we have 
Ana = (1)A, — rA, {27} 


where A, represents the determinant A, with every power sum in the r-th row 
increased by unity. It is only necessary to arrive at some method of designat- 
ing these terms if the above recurrence formula, is to be applied. This can be 
done by inserting the power sum (1) before the other power sums which it is 
to multiply. Also in forming A, add unity to the first power sums in the expan- 
sion of A, being careful to retain the previous order. Thus 


4 = (1) 
& = (0) — 10 +1) = GQ) - @) 
ds = (IGG) — 3) — 2@)@) — By = MMA) — CE) 
— 2(2)(1) + 2(8) 
de = (GG) — (OC) — 20)@)G) + 208) — 3@)0)0) 
+3(2)(2) + 6(3)(1) — 6(4) , {28} 
ds = MMMM — MMME — 2M + 20)0)@) 
~3(1)(2)(1)() + 3(1)(2)(2) + 6(1)(8)(1) ~ 6(1)(4) 
— 4(2)(1)(1) (1) + 4(2)(.)(2) + 8(2)(2)(2) — 8(2)(8) 
+ 12(3)(1)(1) — 12(8)(2) — 24(4)(1) + 24(5) | 


etc. 


By collection of repeated terms and recalling that (1) = A,, the expansion 
becomes 


(1) = (1) 

(1) = (1) — @) 
(1°) = (1)* — 3(2)() + 208) 

(1‘) = (1)* — 6(2)(1)* + 3(2)° + 8(3)(1) — 6(4) {20} 
(1°) = (1)° — 10(2)(1)* + 15(2)(2)(1) + 20(3)(1)(1) — 20(3)(2) 
—30(4)(1) + 24(5). 
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It is possible to write values of (1”) in terms of power sums though the practical 
difficulty increases as r increases. Also continued use of the recursion formula 
{27} is apt to lead to error. Two simple checks are available. If D, represents 
the sum of the coefficients of the expansion of (1’) and | D, | represents the sum 
of the absolute values of these coefficients, then 


D, = 0 when r > 1 {30} 
| D,| =r! {31} 


The proof of {30} and {31} follows directly from {27} since the coefficients 
of A,and A,are the same. Thus D,,, = (1 — r)D, and | D,i:| = (1 + 1r)| D,|. 
Since D, = 0 it follows that D3, Ds, --- , D, = 0 and since | D.| = 2! it fol- 
lows that | Ds; |, | Da |, --- , | D,- | are 3!, 4!, r! respectively. 


26. Determination of the Coefficient of Any Ordered Product of Power Sums 
in the Expansion of A,. We next attempt to revise the process outlined above 
so as to get the formulas {29} without going through the work of writing out 
{28}. We note first that every product of power sums in the expansion of 
(1°) in {28} has been obtained from (1) by a succession of r — 1 operations 
which were either prefixes (when the (1) was prefixed) or raises (when the (1) 
was added). Also the order of the power sums in a given term indicates which 
operations have been prefixes and which raises. For example (1)(1)(1)(1)(1) 
results from 4 prefixes while (5) results from 4 raises. The term (3)(2) results 
from 1 raise, 1 prefix, and 2 raises respectively, while the term (2)(8) results 
from 2 raises, a prefix, and a raise. The product (p.)(p3)(p2)(p:) results from 
prefixes when r = pi, 7 = Pi + p2,7 = Pi + Pe + Ps and raises at all other 
times. 

The sign of the coefficient of (ps)(ps3)(p2)(p1) can be determined when we 
recall that each raise is accompanied by a multiplication by —r while each 
prefix is accompanied by no change in the coefficient. There have been p; — 1 
+p—1l+ps—1+p~—1 = pr + po + Ds + Ms — 4 raises so the sign 
is (—1)" * where p: + p2 + ps + ps =r. More generally if (p,) - -- (ps)(p2)(m) 
is a term in the expansion of (1’) where p, + --- + ps + po + pi = 7 the 
number of changes in the sign isp; -1+m2—1+---+p—-l=p+m 
+---+p—p=r-— op. It follows at once that those products of power 
sums in the expansion of (1’) which have the same number of factors, p, also 
have the same sign and that this sign is (—1)’ °. 

In determining the numerical part of the coefficient we note that each prefix 


is accompanied by a multiplication by unity which can be written in the form -. 


Each raise is accompanied by a multiplication by r-so there appears in the 
numerator the product of all possible values of r and in the denominator the 
product of those values of r corresponding to each prefix. For example the 
numerical coefficient of (p.)(p3)(p2)(p:) is 
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(po + pstmt — 1)! _ (ma+pt+pet+p)!  - 
(ps + D2 + Pi)(p2 + Pi)(pr1) (pa + Ds + D2 + 1)(ps + P2 + Pi)(p2 + D1) (pr) 


Similarly the coefficient, without sign of (p,)(p,-1), --- (ps)(p2)(p:) in the 
expansion of (1°) is 


_— (Pp + Poa + +++ + Ps + m2 + rr)! (32) 
(Dp + Pea + +++ + Pst Pa t+ Pi)(Ppa + °° 
+ ps +p2 + pr) +++ (ps + po + pi)(p2 + Pr) (p1) 
The denominator of {32} has a certain resemblance to a factorial. Thus 
=(14+1+1+41)(01+1-41) (1 + 1) (A) in which the successive factors 
are found by dropping the first unit. The corresponding algebraic expression 


(pe + Ps + D2 + Pi) (Ps + D2 + Pi) (D2 + 21) (p:) is found in the same way 
and might be called an ‘‘algebraic factorial.’”’ It might be designated by 


(ps + Ds + D2 + pridi 


It should be noted that the order of the terms in the algebraic factorial is sig- 
nificant. Thus (pe + p:)i *= (p: + pe)i unless p; = pe. 


The coefficient of (p,)(p,-1) --- (p2)(p:) in the expansion of (1") may now 
be written 


! 
nO aacapaps tenes 33 
ile (Ppp +++ + et Pidi —_ 
For example the coefficient 


(2)(1)(2) is (— Dy —= 6 


(1)(2)(2) is (— "saa 5=3 


(2)(2)(1) is (-* 555 = 8 


and the total coefficient of all terms involving (2)(2)(1) is 15. 
With a less formal notation we might designate the sum of the p! “‘algebraic 
factorials” which can be formed from p, , Pp-1, --+ De, Pi by 


DX (Pp + Peat ++ ++ mdi 
and the sum of their reciprocals by 


1 
Etat + po + pridi 


This notation calls for the inclusion of all the p! algebraic factorials even though 
some of them may be alike. If 


71, My *** 5 Ke 
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indicate the numbers of repeated p’s 


1 
ets ts cats 
1 
se  ciipieninigeineneimienerenin 
ml 2 (Pp + Poi + -°* + 2+ Didi 


’ r! 
where >.’ holds for the na 


!.--+a,! 

In general the total coefficient of (p:)(p2) --- (p,) in the expansion 

of (1’) is obtained by adding all possible terms {33} in which the same p’s 
occur in different positions in the product. Every possible different position 
grouping of the p’s is present but once since it is dependent solely on the unique 
order in which prefixes and raises have been combined to produce that particular 
position grouping. The number of these position groupings varies with the 
number of repeated p’s. The sum of the coefficients of these position groupings 
of the same p’s, i.e. the total coefficient of (p:)(p2) --- (p,) is then given by 


1 
—1)’*r! >.’ —_—_—————— 
(—1)"rl 2 OT Fe 
which can be written by means of {34} 


os r! 1 
(—1) 11! 12! aes eee + pi 


= milas! 


{34} 


non-repeated terms. 









The formula for (1) may be written 


(’) = 2 (-)"” aoa ae SFR (pi)(p2) --- (pp) {35} 


27. Theorem on Algebraic Factorials. 
simplified by the theorem 


1 1 
ee. Sete PpPp-1*** PoP aaa 


which is proved by mathematical indyction. 
A. It is true when p = 2, since 


a ee 
+P (pe + Didi + (pi + Pe)i ste lete 
B. If it is true for p = k, it is true for p = k + 1 since 
1 
—————or 
(Desi + Pe + -°* + pet Didi 











The result {35} can be further 










eet 
Presi tices + Pot pr (pe + ++ + Didi 


1 
° 2» ~ (Pe + DPe-r + °°: isa 
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where 7 gives the k terms in which p+: replaces px, Pet, --* , D2, Pi Te- 
Pkt+1 
spectively. Now if {36} is true when p = k 
1 


a Patatmat tain 
(Pera + °°° 


] 
+ p:)j Prosi tics: + Di | Pk+1 Pk Pe—-1 *** P2Pr 
1 


7 Pe+iPe*** Pop 
C. Hence it is true when k = 2,3,4.---. 


28. Formulas for (i’). Formula {35} may now be written 
oe _ irr? r! 
ay= Dey = i) + @) 
or if the p’s are ordered it may be written as 
r! 1 
") = se FO me ere es... Ee ) 
(1’) a ~ ( 1) my! -+-m,! pit ++ pr (p:) (p.) {38} 
which is the formula previously given as {23} and {26}. In addition the check 
formulas {30} and {31} become 


liv senses =0 {39} 


om! +--+ w,! 


- Ds 
! 
L = >" {40} 


~Pi wi! +++ m,! 


These relations {39} and {40} correspond to statements of Cauchy (2), (I; 1; 
252-3) and to later remarks of Cayley (D; 577). By dividing by r!, they become 


Yq GE = {41} 


p.° m7! eee T,! 
1 
7 1 Ts " = {42}. 
Pi coe De my! ---+ ,! 


The formula {38} is easily — Thus 


. ! ! 2 2 
a) = 26) -[ @@-2O@+ FOO + oy OW 


oe 
— 2 ay’ +3) 
= 24(5) — 30(4)(1) — 20(3)(2) + 20(3)(1)* + 15(2)*(1) 
— 10(2)(1)' + CY 
24 — 30 — 20+ 20+ 15-104 1=0 
24 + 30 + 204+ 204+ 154+ 1041 =5! 


with 
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We next write the formula (1”) in such form that we use the principle of 
generalization from symmetry. If we multiply numerator and denominator 
of {38} by (pi — 1)!""(pe — 1)!" --- (p, — 1)!"* we get 


(1) = © (—1)*(pi — 1) "(p, — 1)!" 


OP DD DD mt a 








which immediately becomes 


(1°) = 2X (-1)*(p: — 1) "(m — 1)! 
’ (ps — 1)!"* La * oe (p:)” re (p.)** {43}. 


This somewhat formidable appearing formula is easy to apply. For example, 
in finding the value of (1°) we write in one row all “a partitions of 5. 
In the next row we place the well known values of * nm a In the next 
row we place the indicated products with proper signs. Thus 
1° 21° 271 31’ 32 41 5 
1 10 15 10 10 5 
1 —1 +1 +2 —2 —6 +24 
results in 
(1°) = (1)° — 10(2)(1)* + 15(2)°1) + 20(3)(1)’ — 20(3)(2) 

— 30(4)(1) + 24(5) 





























as indicated above. 
It is immediately recognized that formula {43} can be obtained from for- 

mula {3} by placing P(1") = (1’); pr* --- pr* by (p:)" --- (p,)"* and Pyn 

by (—1)" °(p,1 — 1)!"" --- (p. — 1)!"* and hence that formulas of Table I 

may be used in obtaining the values of (1’). 












29. Values of (a; --- a-). The form of {43} also permits generalization 


from symmetry since the * , pt equal values (p:)"' --- (p.)"* are re- 
— 


& 
r 


placed by the Lia . pt different values composing T(p,)"! --- (p.)"* when 


the r units are replaced by the a’s. It follows at once that 

(ayaz --- a,) = 2(— 1) *(p, — 1)!" --- (pe — 1) IT (p:)™ --- (p.)** {44} 
where r = Pim, + Powe +--+ + Pete 
and p m+ mMet+--- + mm. 






















li 
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As an illustration we write 


(abc) = 2(a + b + c) — (a + b)(c) — (a + c)(b) — ( + c)() 
+ (a)(b)(c) 


as indicated earlier by {18} and 


(a,a2a;a,) = — 67T(4) + 27(3)(1) + 7T(2)(2) — T(2)0Y¥ + TA) 
(aya2q30,a5) = 247(5) — 6T(4)(1) — 27(3)(2) + 27(3)(1)’ 
+ T(2)’Q1) — T(2)Q)* + TA)’ 


ete. 


30. Table of Values of (a; --- a,). The values of the power products with 
w =< 6 are given in Table II which follows the general form of Table I. In 
fact Table II may be derived from Table I by placing every 


Pyn...9%8 = (— 1)" “(m — 1)!" --- (p, — 1)!" 
as indicated in the next section. 


31. Use of Partition Formulas. By comparing {44} with {4} we see that 
{44} can be obtained from {4} by placing 


P(a,a2 ve a3) = (a,a2 aia a,) 
Pyq...pte = (— 1)" *%(p. — 1)!" --- (pe — 11" 
and T yq1...pte = T(pi)™ --+ (p.)" 


It appears then that the values of any power product sum (a;a2 --- a,) can 
be obtained by writing the expansion of P(a, --- a,) and substituting as indi- 
cated. Thus since 


P(321) = P36 + Po51 + P42 + Py33 + Py 11] 
(321) = 2(6) — (5)(1) — (4)(2) — ()@) + ()G)Q). 


It is also immediately apparent that Table II can be obtained from Table I 
by placing P,7:...,r. equal to (— 1)" *(p,: — 1)!" --- (p, — 1)!"* and that the 
main results of Chapter I, including the recursion rule, are applicable to the 
present problem. 


32. Coefficients of Given Terms in the Expansion of Product Power Sums. 
The methods of the last section are also useful in finding the coefficient of any 
term in the expansion. For example we wish to find the coefficient of (3)(2) 


in the expansion of (2111). We note that a) = Px, + 3P2x2 and that 
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TABLE II 


Power product sums in terms of products of power sums when W S 6 
W=6 
we 
33} 411 | 321) 222)3 2212 |214 | 16 



































| 45 |-15 


5| 41| 32 (| 311 2111|11111 





























the coefficient of (3)(2) is Ps: + 3P22 where Ps, = (— 1)*°2! = 2 and Py 
= (— 1)** =1. Hence the coefficient is 1-2 + 3-1 = 5. 


33. The Expansion of the Monomial Symmetric Function. If a; a2. as; 
% --- a, then M(q --- a,) = (a --- a,) and previous results are applicable. 
If however the product power sum is of the form 


(af! az? on ay") 
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1 


M(az" az? --- ax") = —— 
@1.@q. °° 


, (ar* +++ ay”) 
* Ae 
and 


M(az' a% 


aj") 





: a (—1)" "(pr — 1)" --- (ps — 1)"* TU)" --- (p.)"* {45} 


a! axe! 
For example 
M(421) = 2(7) — (6)Q) — ()Q@) — 4G) + @@)Q) 
M (322) = (7) — (5)(2) — 2(4)G) + 2@)@)(). 
M(2'1*) = —3(6) + (5)(1) + $(4)@2) + 3)@) — 44)()Q) 
— 3)@2)Q) — #2)@)Q@) + 22)2)G)Q). 


Study will show that the formula {45} is equivalent to one given by Faa de 
Bruno (C; 9) and later by Roe (7). 

It is possible to.use Table II in finding the expansion of the monomial sym- 
metric functions. It is only necessary to multiply each term in the expansion 
of (a; --- a,) by sesgettmnnan, 

Gai °°* Gat 
The check formulas give, in the case of the monomial symmetric function: 
The sum of the coefficients in the expansion is 0. 


! 
The sum of the absolute values of the coefficients is ——”~——_.. 
; a! a! 2 oe a! 

The reader might compare the second of these checks with the results of 
Faa de Bruno (C; 14). 


Tables giving the expansion of monomial symmetric function have been 
given. One by J. R. Roe (12; plate 18) includes all cases of weight <10. 





34. Previous Results. Previous authors have studied the monomial sym- 
metric function. Gordan has deduced a monomial symmetric function formula 
which is recommended by J. R. Roe (M; 24-33). MacMahon has given a 
general formula (K; II; 320) for expanding any monomial symmetric function 
in terms of power sums together with an operational method for its evaluation. 
O’Toole also has given a differential operator and showed how it could be applied 
in obtaining expansions (16; 115-130). O’Toole has also given a method of 
expanding symmetric functions in many variables by means of differential 
operators, (17). 

Another method of attack was based upon the close relation existing between 
the elementary symmetric function and the determinant of the power sums. 
This has resulted in the expression of the monomial symmetric function in 
determinant form. Brioschi appears to have been the first (1854) to see how 
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a symbolic determinant could be used (3; 427) although he gave no proof. 
Bellavites tried in 1857, but obtained incorrect results (4). In 1876 Faa de 
Bruno made an attempt, but he too was in error (C;10). In 1898 E. D. Roe, Jr. 
proved that Brioschi was right (7). Muir also gave a proof in 1908 (11; 5-9). 
The summation of determinants, rather than the symbolic determinant, was 
used by Hankel (6; 90-94) (L; III, 220). 

The determinant of the power sums has been generalized in another way. 
A group of writers has studied the “immanents” of its matrix. D. E. Little- 
wood and A. R. Richardson have recently written a series of papers on this 
topic. One of these papers (18; 99-141) defined the term “immanents” and 
gave references to previous investigations dealing with this matrix. 

It has been the aim of this chapter to present an easy development of the 
subject of the expansion of product power sums and monomial symmetric func- 
tions. This development is characterized by 

1. The use of the formulas and tables of Chapter I in writing expansions of 

product power sums. 

2. The use of product power sums in place of monomial symmetric functions 

which makes feasible 

3. Generalization from symmetry. 

4, References to previous work. 


Chapter IV. The Double Expansion Theorem 


In the present chapter we combine the multiplication expansion of Chapter II 
and the power product sum expansion of Chapter III into a new result which 
is to be known as the double expansion theorem. We show that this result 
may also be expressed in terms of the partition notation of Chapter I. 

¢ 


35. The Value of K(a;)(a2). We know 


(a:)(a2) = (a; + a2) + (ai-a2) 


and if we multiply (a; + a2) by ke and (a;-a2) by ky, we have a new expression 
which we designate by K(a;)(az). 


K(a:)(a) = ke(ar + a) + ku (ar-a2) {46} 
Since (aiaz) = (a:)(a2).— (ai + ao) 

K(a:)(a2) = ke(ai + a2) + kul(a)(@) — Ga + a)] 

K(a:)(a2) = (ke — ku)(ai + a2) + ku(ai)(@) 
which can be written 

K(a:)(a@2) = Ke(ai + a2) + Ku(ai)(a) {47} 
if ke — ky = Keand Ky = ky 


f 


iS 


I 


It 


on 


17} 
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36. The Value of K(a;)(a2)(a;). We know from {12} that 
(a:)(a2)(as) = T(3) + T(21) + T(11) 


and we define K(a;)(a2)(as) = ksT(3) + ka T(21) + kinT(111). Inserting 
the values T's = (a; + a2 + as), Tn = (a, + ae-a3) + (a, + 3-2) + (a2 + a3-a), 
Ti = (a,a2a3) and reducing to power sums by {44}, we get 


K(a:)(a2)(as) = (ks — 3k + 2kin)(ai + ae + as) + (kor — kin) 
{ (a: + a2)(as) + (a: + as)(a2) + (a2 + as)(a1)} + Kin (ai)(a)(as) 
which may be written 
K(a;)(a2)(as) = Ks(ai + a2 + a3) + Ku{(ai + a2)(as) 
+ (a: + a3)(a2) + (ae + a3)(a1)} + Kin(ai)(@)(as) {48} 
where K; = ks — 3kq, + 2kin, Kn = kn — kim, Kin = kim. 








37. Definition of K(a;)(a2) --- (a,). We define | 
K(a;)(a) --- (a,) = a kpth...pte T(pi' --- p:') {49} 


where T(pi' --- pz*) is composed of { _, : x, ) power product sums. We 
pi D sims Ds 


wish to find the value K(a;)(a2.) --- (a-) in terms of power sums. This in- 
volves the expansion of each power product sum in terms of power sums and 
then the collection of the results. This algebraic process is to be called the 


double expansion process and the theorem which results, the double expansion 
theorem. 


38. Special Cases of the Theorem. The results {47} and {48} are special 
cases of the double expansion theorem when r = 2,3. Whenr = 1 it is evident 
that K(a:) = Ki(a:) = h(a). {50} 
The results {50}, {48}, and {49} may be written symbolically by 


K(a) = KiT(1) 
K(a;)(a2)(as) = K;T(3) + KuT(2)Q) + KwT(1)’ 


It can also be shown, with a much more extensive use of the results of Chapters II 
and III, that 


K(a:)(a2)(as)(a) = K.T(4) + KuT(3)(Q1) + KuT(2)’ 

+ KoyT(2)Q0) + KunTAy {52} 
K(a) --- (as) = KsT(5) + KaT(4)Q) + KnT(3)(2) 

+ KsuT(3)(1)’ + KeuT(2)°Q) + KemT(2)(1)’ 

+ KumT(1)° {53} 
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Ky = ky — 4k, — 3k + 12ken — 6hin 
Kau = ks —3hkou + 2kun 
kee — Qkeu + kim 













Kou re kun 


kun 


Ks = ks — 5ka — 10ks2 + 20ksn + 30ke21 — 60k + 24k 
Ky = ky — 4ks. — 3k + 12kean — 6kinn 
Kz = ks — =k — Skea + Skoan — 2kian 
Ksu = Kau — Skeaun+ 2k > {55} 
Ko, = Kooy — 2kernr + 
Kon = 
Kian = 












































Kun 








kon WT: Kun 










kinu 


We may say then that, for r < 6 
K(a:) --- (ae) = Do kop..-pte T(pi! ++ pe) 
= Do kpp...ote T(ps)" «++ (p.)"* {56} 


where K,r1...,%4s defined by the relations {47}, {48}, {54}, and {55}. In 
examining the value of K, we note 


Ki =k, 
Ke = ke — ky 
K; = ks — 3kn + 2kin 
Ky = kg — 4ka, — 3ko + 12k — 6hin 

Ks = ks — 5ka — 10ks2 + 20ksn + 30K. — 60Kan + 24h 
and that these are given, for r < 6 by 





















+ 


pe — i Tl Ls 
K, = 0 (-1)(0 D!( wel Keptt...pte {87} 


It is further to be noted that {57} can be obtained from {3} by placing P(1’) 
= K,, pi’ --- ps’ by kyn...pte, and Pyn...p%8 by (—1)’(p — 1)! Hence the 
last rows of Table I may be used in writing the values of K,. Thus from 


P(1’) = P;(3) aod 3P2(21) + Py(1)° 
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we get 

Kz; = ks — 3ke + 2kin. 
It is further evident that if K3Ke = (ks — 3ke + 2kin)(ke — ky) indicates 
multiplication by suffixing of subscripts that K3;K2 = ks. — ksu — 3kean + 5keu 


— 2k = K3e 
and in general it can be shown that for r < 6 


Kr, = Kr, Ke, 
Kyyrr3 = Kr, Kr, Ke, {58} 


ete. 


so that all values K,1...,%. may be obtained by symbolic multiplication of 
equations {57}. 

The method of this section can be used in demonstrating that the results 
{56}, {57}, and {58} hold also when r = 6, 7, 8 --- , but the amount of alge- 
braic manipulation increases enormously with each increase in r. We establish 
these results, for all integral values of r, by a more general approach. 


39. A More General Definition. We provide a more general definition of 
K(a;) --- (a,) by letting the subscripts of the k’s agree with parts of the given 
partition rather than with its complete order. Thus 


K'(a;)(a2) = Kea, +0,(a1 + 2) + Kaya,(a102) 


and in general, if gi’ --- gi‘ represents any p part partition having complete 
order pj' --- p,*, then we may define 


K'(a:)(az) «++ (ar) = Do heogt..-azt (qi' - ++ ae*) {59} 


where the summation holds, not only for every different complete order as 
does {49}, but for every possible partition. By {44} (qi' --- gi!) may be 
written as 


Gi --- gi) = D(-1)" “dG — )M@ — 1)! --- @ — IT) --- @) 


where d; + dz + --- +d, = p and where groups of the d’s may be alike. If 
(w:)(w2) --- (w,) is one of the products of power sums having the complete 
order (d, --- d,) we may write 


(qi! --- a) = Do(—1)" “(di — 1)! --- G@ — 1) (ws)(w») --- Cw) ~——( 60} 
where 


M21 + --- +t = W= UW + We +--+ + Uy 
and 


m+:--+2%=p 
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and where the summation sign holds not only for every complete order d, - - - d, 
but for all power sum partition products (w:)(we) --- (w,). 
The insertion of {60} in {59} gives 


K'(a;)(az) - - - (a) 
= Zz Keg. .gZt Z (—1)’ "(ad — 1)! vee (d, = 1) !(w) = (w,) {61} 





40. Value of K.,. The notation of K, is used to indicate the coefficient of . 
the power sum (w) = (a; + a2 + --- + a,) in the expansion of {61}. In this 
case d; = p and g = 1 so that 

Ki, =  kegzt...qze (—1)" “(p — 1)! 
which may be written more symbolically as 
Ky =  (-1)" "(0 — I)! ks, {63} 


where 7, represents any algebraic partition of a; + --- + a, and p indicates 
the number of its parts. 


{62} 





41. Products of K’’s. The notation Ki, Kos is used to indicate the product 
of Ki,, by Ki, if the rule of multiplication is the suffixing of the subscripts of 


the k’s in the expansion of x... and ol . Thus 
niin — | — i 
= Kay+ay-a3 me Kayazay 


More generally, if we write, from {63} 
Ke, = 2 (-1)* "(hh — I! he, 
Ko. = Do (—1)7 (de — 1)! eee, 


SCHOKSHSE SHA SESE HCH SECCHEDR BHCC O ST HHS 


Ki, = (1), — 1)! kee, 
and use multiplication by suffixing of subscripts we have 
Ko, Koy +++ Kg = 25 (—1)" (di — 1)! - + (dy — 1) beng, ---0, 


where p = dg, + d& + --- + d, and the summation holds for every partition 
which can be formed by combining any algebraic partition of w, , any partition 
of we, --- , any partition of w, . 


{64} 


42. The Coefficient of (w:)we) --- (w,). The coefficients of any specific 
product of power sums (w;)(we) --- (w,) is from {61} 


Kooywy---0, = >. keggt...gze (—1)° (di — 1)! --- (dy — 1)! 





{65} 
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where the summation holds, not only for the partitions of a; + a, + --- + ,, 
but for the partitions tu,, Tw,, ---, Tw, Since these partitions can be combined 
to form (w;)(we) --- (w,). Henee {65} becomes 


/ 


Ky wg: ++, _ Zz (—1)’ "(a - 1)! a (d, = 1)! Kiry,++-t0, {66} 
and it is immediately seen that the right hand expressions of {66} and {64} 
are the same and hence that 


, _* ——--—~”~!lhlUuuee 
iain Bale + 


g 


as expected from (58). 
We can now say that 


K'(ai)(a2) «+ - (ar) = DO haz. .age Qi! = ++ gf) 
= DY Kujwe---w, (wi)(we) - + (w,) {67} 
where 
Ky = D (-1)""(o — It ke {68} 
and 
Roig * Balen *** ee {69} 


Relations {67}, {68} and {69} constitute the general double expansion 
theorem. 


43. The Double Expansion Theorem. The case of the double expansion 
theorem in which we are especially interested is that in which the coefficients 
of all similar power sum products are the same, i.e., k,z1...,z¢ is a function of 
the complete order indicated by kyt...,%+. In this case {68} becomes 


— _— ° — 1 1 | T1 , 
Ke =X (1 = (pos es) bates (70) 


where the summation holds for all possible complete orders. Suppose now 
that the r algebraic expressions, a , d2, --- , ad, are all unity then {69} becomes 


” 1" 
K, = Z (2%, D1 wh on) Kpr...pte 
and we find that K. = K,. We may then write {67}, {68} and {69} as 
K(a,) --- (a,) = } Kpit...pse (pi' --- p:’) - zz Ky,...1 T(r) “0% C) {71} 


where 


K= D(-v- (nT (72) 


iii Ma {73} 
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Now rrz --- 7, indicates any grouping of the a’s, and hence any complete 
order of a; + a2 + --- + a,. So {71} may be written, with a slight change 
of notation as 


K(a;) oe (a,) = Zz. Kpt1...pte T(pi' — p:') 


= 7 Kar...pre T(p:)" ve (p.)”* {74} 
The relations {74}, {72} and {73} are the desired generalizations of {56}, {57} 
and {58} and hold for all positive integral values of r. 

The double expansion theorem provides a method of writing out the result 
of the double expansion process without going through the work involved in 
the process. Thus 

K(3)(2)Q) = K3(6) + Kn{(5)G) + ()@) + @)G@)} 
+ Kin(3)(2)Q) = (ks — Bho + 2kin)(6) {75} 
+ (ka — kin) {(5)Q) + 4)Q2) + B)GB)} + kn B)2)Q) 





44. The Double Expansion Theorem and Partition Notation. It is im- 
mediately evident that {74} can be obtained from {4} if P(a, --- a,) is replaced 
by K(a) --- (a,), if Ppr1...,% is replaced by K,n...p%, and if pr’ --- py*is 
replaced by T(p:)"' --- (p.)"*. It follows at once that the entire theory of 
Chapter I,—table, recursion formula, etc.—is applicable to double expansion 
theory. For example {75} above is obtained from 


P(321) = P36 + Pa{51 + 42 + 33} + Pin321 
simply by replacing the K’s by the P’s and enclosing the parts in parentheses. 
We can as well use P’s as K’s to represent the double expansion theorem and 
hence have available a list of double expansion formulas when w S 6. We 
also have available a recursion property for writing double expansions beyond 
the scope of the table. Thus for éxample, the illustration at the end of sec- 


tion 9 may be interpreted as a statement of the double expansion theorem 
when a; = 3, a2 = 2, a3 = 2, a, = 1. 








45. The Case of Equal Powers. In case a; = a2 = a3 = --- = a, {74} 
reduces to {3’} of Chapter I with 
a 1’ 
P, = —1)*" -1(,, +.) katteot 76 
Zz ) (p ) Pi'---p,° Pi Ps { } 
and 


Poiaicig Pin «+ Py. 


Formula {74} also reduces to {3} when a; = a, = --- 






46. Special Values of K,71...,%+. 
A. kpn...pre = 1. In this case the coefficients are all unity and 


P(a)(az) --- (ar) = (a)(@) --- G) 
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It follows that P, = 0 and that P,1...,7. = 0 except that P, = 1. Placing 
P, = 0 and k,n...,r. = 1 in {72} or its equivalent {76} we have, when r > 1 


0= Dy w- I(T Le.) (7} 


where the summation holds for every partition of r. This formula should be 
compared with {39} and {40}. When, = 4 and the partitions are 


4, 31, 22, 211, 1‘ 
{77} gives 1 —4 —3 +12 —6 
{39} gives -—6 +8 +3 —6 +1 0 
{40} gives 6 +8 +3 +6 +1 = 4! 
The equivalent of {77} was first given by Cayley (D; —_ who at the same time 
noted the similarity to {39}. 

It follows immediately that the sum of the coefficients in the expansion of 
P,11...pt*, except Pir, is 0, for the sum of the coefficients of P,z1...,%. is the sum 
of the coefficients of (P,,)" --- (P,)"* and is 0. For example the sum of the 
coefficient of P32 = kee _ ku _ Skee om 5ke111 _ Qh is 0. 

Since the coefficients of (ui)"! --- (u.)"* (19; 25) in the expansion of Thiele 

- 


half invariants are (—1)”‘(p — 1)! Ln ol it follows from {77} that the 
Ps. 6a 


sum of these coefficients is 0. 
(p) 


B. kpz...p72 = a In this case all terms having the same number of parts, 


n” 
p, have the same coefficients. If we indicate —— by pi, pe, ---, whenp = 1, 


N® 
- , {57}, {76} become 
Py = pr 
Pe = pi p2 
Ps = pi — 3p2 + 2p; 
Py = pi — Tp2 + 12p3 — 6p 
Ps = pi 15p2 + 50p3 — 60p, + 24p;. 
etc. 


which are the formulas which have been used by Carver (15) and O’Toole (16). 

Many other additional cases can be obtained by giving different values to 

nl ..pt#, but a discussion of these is hardly justified here as the case in which 
.pts is a function of the number of parts, p, is to be used in Part II. 


47. Relation to Previous Results. No general statement of the double 
expansion theorem has previously been given although the special case K(a’) 
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has been developed by Carver (15) and O’Toole (16). Their results are further 
restricted to the special case (B) of section 46. The application of the double 
expansion theorem in this case is very useful in studying sampling from a 
finite universe as Carver has shown and as is demonstrated in Part II. 

Most writers who have worked on the problem of moments of moments have 
gone through the double expansion process, but Carver was the first to note 
that the result of the process can be written in terms of the P polynomials 
above. It seems appropiate therefore to refer to these P polynomials of the 
coefficients as Carver polynomials. 


Chapter V. The Multipartition and Multivariate Formulas 


It is the purpose of this chapter to show how the results of Chapters I, II’ 
III, and IV may be extended to the case of different variables. 


48. Multipartitions. ‘ables. Formula {4} is still applicable if we let the 
a, units be the units of one quantity, the a2 units to be the units of a second 
quantity, etc. Thus for example the formula P(a,a,a3) may be used to repre- 
sent the precise number of ways in which a; apples, az pears and a3 peaches can 
be formed into groups without breaking up the groups of apples, pears, and 
peaches. 

Various conventions for representing multipartitions of this type have been 
used. We adopt the one in which the individual partitions are written in 
successive columns. The partitions of the first number are combined with the 
partitions of the second number to form all possible multipartitions. Thus the 
multipartite number 111 has the partitions 


111 110 101 011 100 
001 010 100 010 
001 


where the parts are given in the rows. It is desired to show the number of 
ways in which any-one of those partitions may be combined to form partitions 
of fewer parts. Thus 


100 110 101 011 100 
010 } = P3111 + P2001 + P2010 + P2100 + P,010 


001 001 


This is obtained from P(a,a2a3) by placing a; = 11, dg = 12, as = 13, and could 
be written from {4} as 


P(1i1213) = Ps(1i + 12 + 1s) + Poifdy + le-13 + 11 + 13-12 + le + 13°11} 
+ Piyli-1s-1s. 
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Similarly 
10 20 02 
P 10 = P,2 aa 2P3,21 + 2P3,12 a P2220 ae 2Py411 oe P 2,01 a P2o,10 
01 2 01 10 02 11 01 10 
01 
11 10 
+ 4Po10 + Punl0 
01 01 
01 


is a special case of P(a,a2a3a4) where a; = 11, d2 = 11, us = le, a, = 1g. For- 
mula {4} is also true where the a; units are not of the same kind. Thus 


P(a,a2) = Pe(a,; + a2) + Piu(aiaz) 
gives 


P(io) = P221 a Puto when a= 1, + le and a= 1. 


TABLE III 
The Multipartite Number 11 


| 10 
01 





Pu 





The Multipartite Number 111 


11 | 110 | 101 | O11 
| 001 | O10 | 
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TABLE I]I—Continued 
The Multipartite Number 22 


12 | 20 | 11 | 20 | 02 
10 | 02} 11 | 01 | 10 

















2Px | 2Pu | Px | 2Px Pm | Pau | 
Poy | 


| 
| 














Tables can be made for the partitions of the various multipartite numbers. 
In Table III are presented values for the numbers 11, 111, 22. 
When the units are indistinguishable 11 condenses to the w = 2 part of 
Table I. 
When the units are indistinguishable 111 condenses to the w = 3 part of 
Table I. 
When the units are alike 22 condenses to the w = 4 part of Table I. 


49. Multivariate Distributions. The chief results of Chapters II, III, IV 
also hold for multivariate distributions. Some additional definitions are neces- 
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sary. We suppose that the N variates z,, a2, --- , 2y are replaced by the 
Nr variates of the array 


1%1 , 142, °** »1%N 


o%1 , 2%2,°°° » 2UN 


rX1, rX2, ale » rUN 


where the presubscript represents the variable. The power sums become 
a 
(a;) = yy" + 123" oe) + av = >» 1x5" 
(a2) = oxi? + or3? + --- + or? = Do ox? 


It is not necessary to utilize the presubscript since it is precisely the subscript 

of the a. That is the power sum (a,) is defined by >> z{*. Similarly (a, a2) 

= >> 123" ex? can be written as (a;a2) = >, 2%! x3? without introducing ambiguity. 
oj tej 


In general {6} as well as {4}, now holds for the multivariate case. It follows 
at once that the results of Chapters II, III, IV can be written for the multi- 
variate case by means of the formulas of Chapter I as indicated by the previous 
section. Thus the formula for P(1,1,l212) may be written as [Table IIT] 


P[10-10-01-01] = P,22 + 2P21-01 + 2Px112-10 + P22002 + 2P 1111 
+ Pon 20 01 01 + P2102 10 10 + 4Pen11 10 01 + Puni0 10 01 01 
and can be interpreted as: 
(10)°(01)° = (22) + 2(21-01) + 2(12-10) + (20-02) + 2(11-11) + (20-01-01) 
+ (02-10-10) + (11-10-01) + (10-10-01 -01) 
by {12} of Chapter II. It may also be interpreted as 
(10-10-01-01) = — 6(22) + 4(21)(01) + 4(12)(10) + (20)(02) 
+ 2(11)(11) — (20)(01)(01) — (02)(10)(10) 
— 4(11)(10)(01) + (10)(10)(01)(01) 
by {44} of Chapter II. It can also be interpreted as a double expansion by 


means of section 44 where the values of the P’s are given by the usual 


= 1’ 
PF. => —] — 1 ! Fi v1) Batt Te 
EM —— (es Pes) hata 


ee 


Pyy-org = PyyPry +++ P 


T9 


50. Summary. It is apparent that {4} not only expresses (a) the number 
of ways in which the parts of one partition may ve collected to form the parts 
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of another partition, (b) the formula for expanding products of power sums 
in terms of power product sums, (c) the formula expanding power product 
sums in terms of power sums, and (d) the formula for double expansions, but 
also that it can be used to make similar expansions in the case of multivariate 
distributions. 
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ON THE INDEPENDENCE OF CERTAIN ESTIMATES OF VARIANCE’ 


By Atuen T. Craic 








1. Introduction. It is well known that a necessary and sufficient condition 
that several statistics be independent in the probability sense, is that the char- 
acteristic function of the joint distribution of these statistics shall equal identi- 
cally the product of the characteristic functions of the distributions of the 








individual statistics. Thus, if 27, 22,---,2y are N independently observed 
values of a variable xz which is subject to the distribution function f(z), and if 
6, , 02, --+-, 0, are s statistics, each computed from the N observed values of z, 






the characteristic function of the joint distribution of the s statistics is given by 


(th, te, _ zs - / ae | ee -omnge) Deca f(rw) dxy Sahn dz. 


Here, i = ~/—1 and the limits of integration are taken so as to include all 
admissible values of zx. Since the characteristic function of the distribution of 
6,,v = 1, 2, --- , s, is given by 


aie / ee ‘ eo f(x.) ... flay) day -» dan, 


the necessary and sufficient condition for the independence of the s statistics 
can be written 


(1) o(ti, +++, te) = gilt) --- golts), 
for all real values of t, , 2, --- ts. 


An important phase of sampling theory in statistics is that in which the 
variable z is subject to the normal distribution function 
22 


1 wie 
ee ee —-eSzt<¢-%, 










oV 2r 












and 6,,---, 6, are s real symmetric quadratic forms in the N independently 
observed values of x. That is, 





Me 
Maz 


~ 
li 

~ 

> 
I 


6; = Aj, Lj XK , 












1 


Mz 
Mz 


= ( bine UjXe , 


~ 
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— 
> 
H 
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6, = 





Ms 


Pik Vj; Xk , 
1 






1 Presented to the Institute of Mathematical Statistics on December 30, 1937, at the 
invitation of the progam committee. In the paper, we discuss, from a slightly different 
point of view, some of the material found in the references given at the close of the paper. 
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so that 


ee ee ee 


where T = t zy ke Gjntjte + --- +t, Z be Dik UjX~ — 5 i. xj. If 


Ai, +--+, As denote the real symmetric matrices of the s quadratic forms, the 
characteristic function can be written 


o(ti, --+,t) =| I — 2iot,Ai — .-- — 2io*t, A, |, 


where I is the unit matrix of order N and the vertical bars indicate the deter- 
minant of the matrix within them. Similarly, the characteristic function of 
the distribution of 6, is given by 


gts) = | I — 2io*t, A, |, 


so that a necessary and sufficient condition for the independence of the s real 
symmetric quadratic forms can be written 


(3) | I — 2io?t,Ai — --- — 2io’t, A,| = [| 1 — 2i0°t, Ao |, 
v=1 


for all real values of t) , 2, --- , ts. 

Although equation (3) is fundamental and is of considerable value in certain 
problems, it should be remarked that it is frequently rather tedious to use. 
This suggests that by strengthening the hypotheses, it may be possible to 
establish another necessary and sufficient condition which, in certain cases, 
may be easier to use. 


2. Certain quadratic forms. In order to lead up to such a theorem as that 
suggested at the close of the last section, we first consider two theorems regard- 
ing real symmetric matrices. . 

Theorem I. ‘Let Ai, Az, ---, As be s real symmetric matrices, each of order N, 
such that Ay + Ao + --- + A, = I, where I is the unit matriz of order N. Let 
tr ,v = 1,2, --- , 8, be respectively the ranks of the matrices A,. If: +72+ --- 
+r, = N, each of the non-zero roots of the characteristic equations’ of the matrices 
A, is +1. 

If s = 2, the theorem is almost self-evident. For the characteristic equation 
of Az is | A2 — AJ | = 0, which, since A; + Az = I, can be written | J — Ai — 


2 By the characteristic equation of the square matrix A is meant the algebraic equation 
of degree N in X.|A — AJ|=0. If A is real and symmetric and the rank of A is r, the 
characteristic equation has exactly r real non-zero roots and N — r zeroroots. Cf. Kowa- 
lewski, Einfihrung in die Determinanten-Theorie (1909) pp. 126-128. 
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I | = Oor| A; — (1 — A)I | = 0. But the last equation is the characteristic 
equation of A, with \ replaced by 1 — X. Thus the roots-of the equation 
| A, — AJ | = O are one minus the roots of | Az — AJ | = 0. Since the equation 
| As — AI | = Ohas N — 72 zero roots, the equation | A; — AJ | = OhasN — ry, 
roots equal to +1. But 7, = N — r so that all the non-zero roots of | A; — 
J | = Oare +1. A similar statement holds for the roots of | A2 — AJ | = 0. 

In general, we have A, + A2o+---+A,=Jandn+n+---+7=N. 
Let B, = Ae + As + .--- + A, and denote by R, the rank of B,. Thus’ 
Ri <re+7r3+--- +7,. Now A, + B, = I and the equation | A; — J | = 0 
has exactly N — r,; zero roots. Since the roots of | B, — J | = 0 are one minus 
the roots of | A; — AJ | = 0, the first of these two equations has at least N — 1, 
non-zero roots so that Ry > N —m=re+r3+--- +7. Fromr+73+4 --- 
+r,< Ri S re +73 + --- +7, we deduce the equality so that the argument 
in the case of s = 2 applies to the matrices A; and B,. In particular, then, 
each of the non-zero roots of | A, — J | = 0 is +1. By writing B, = A; + 
A3;+---+A,, Bs = Ai + Ao + As + --- + A,, and so on, and repeating 
the argument in each instance, we see that the theorem holds. 


Theorem Il. Let Ai, A2,--+, As be s real symmetric matrices which satisfy 
the conditions of Theorem 1. There then exist s — 1 real orthogonal matrices of 
order N, say L, , Le, --- , Ls. , such that each of the s matrices 

Rg ++ BEAM, ooo Bee, v=1,2,---,8, 


is a diagonal matrix‘ with the r, non-zero elements on the principal diagonal equal 
to+1. Necessarily, the sum of these s matrices is the identity matriz. 

In proof of the theorem we shall, to save space, restrict ourselves to the case 
of s = 3, although the method we use will be readily seen to be entirely general. 
Since A, is real and symmetric and since, by Theorem I, the r; non-zero roots of 
the characteristic equation of A; are +1, there exists a real orthogonal matrix of 
order N, say L; , such that 
Bice. 
eer : 


1 : 
0 0: 
















where L; is the conjugate of L,; and where, merely as a convenience of notation, 
we have placed the 7; non-vanishing elements of the principal diagonal in the 
first 7, rows and columns. If then, in both members of the equation A; + Az + 
A; = I, we multiply on the left by L; and on the right by L; , we have 






3 Cf. Bécher, Introduction to Higher Algebra (1921) p. 62. 
‘ By a diagonal matrix we mean a matrix whose elements not on the principal diagonal 
are zero. 
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- 0 


+ LiA2L, + Li Aska = I, 


since Li IL, = IL;L,; = I. The matrices L}A2L; and Li A3L; are real, sym- 
metric, and the ranks are re and r3, since L; is non-singular. Moreover, the 
non-zero roots of the characteristic equations of the two matrices are +1; for 
|LiAeI, — | =|Li(A2 — A)L,| =| Li || A2 — AZ || Li|, and similarly 
for the matrix L;A;L,. Now if a real symmetric matrix is positive definite, 
that is, if all the non-zero roots of its,characteristic equation are positive, then’ 
all the elements on the principal diagonal are positive or zero, and, if an element 
on the principal diagonal is zero, all the elements in the row and column in which 
that element lies are zero. These two facts regarding a real symmetric positive 
definite matrix, in conjunction with equation (4), require that the matrices 
L;AzL, and L;A3l, be of the forms 


» 6 bri4tri4l ot* con® Crittrptl * °° Cry4iN 


- 0 bwiry4t soe @ Cyyri41 **+* CNN 


respectively. Now the real symmetric matrix 


be 4tirigt °° Or 4a. 
C= . 
bwyry4t 
is of order N — r,, its rank is re, and its characteristic equation has r, roots 
equalto +1. There then exists a real orthogonal matrix M of order N — 7 , say 


Mery 4tyry4l °** Mery 41,N | 
1 


M 


such that 


5 Cf. Cullis, Matrices and Determinoids (1918) vol. 2, p. 302. 








52 ALLEN T. CRAIG 


Again, to simplify the notation, we have placed the 72 non-vanishing elements of 
the principal diagonal in the first rz rows and columns. Consider the orthogonal 
matrix of order N 


1 0O - 0: ere ee 0 
0 1 - 0: 
Ig = 0 O--- 1) O-- seen eee 0 
0 i ees ee 0 Mr 41,7141 Mr +1 N 
Q.----- QO Myyry41 ee: ™nNn 


It is evident that L(L1AiL1)L_ = Li; AiL;. If then, both members of L;AiL; + 
Ad. + L, A3L, = I are multiplied on the left by Lz and on the right by L., 


we get 
0---00------ @...6 
1 0---00---0 i 
i 0--- 0: 0---0 
ore <a) oe 
sph ahsicdaindleeel + ]0---0:1 0 0 
rere 00 0 vie 
eens 
A) sbsoneeewc sen 0 : 
0 icccnmaiamews 0 
0---0: O---eeeee 0 
0 0 ree 0 
Becta =] 


0 ssi 0 dy s44,7341 ae te 


Q0-.--0 dy,r,41) -+++ dyy 


From this last equation, it follows that dj, = 0,j7 # k,d;; = 0,7 =n +1, 
-2+,% + re and d;; = 1,j =m +7+1,---,N. The third matrix in the 
left member of preceding equation then takes the form 


0 --- 0: O-seeeeeeeeeee 0 
0.-.--0 
0 0---0 
0-..0. 0 
1 0O.---0 








of 
al 


~~ 
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This establishes Theorem II when s = 3. The procedure may be continued 
in a fairly obvious manner so as to justify the theorem for any finite positive 
integer s. 

With the aid of Theorems I and II, we are now able to state and prove a very 
useful theorem on the independence of certain quadratic forms of normally 
and independently distributed variables. The theorem follows. 

Theorem III. Let 21, t2,---, ty be N independent values of a normally 
distributed variable x and let 0,, A 6, be s real symmetric quadratic forms in 


these N variables, where Zz 6; = Zz aj. If r,, re, +++, Ts denote respectively the 
1 1 


ranks of the quadratic forms, a necessary and sufficient condition that the s forms be 
independent in the probability sense is thatr; + re + --- +7, = N. 

Consider the characteristic function of the joint distribution of the s forms as 
given by equation (2). In accordance with Theorem II, we can successively 
introduce new variables by performing real linear transformations with orthog- 
onal matrices L,, L2, --- , Ls; respectively in such a way that® 7 becomes 


Tr} ritre N 1 N 
T=hDyith Dyit--- +h Vi — 5 Die 

I rifl rte brit 20? “F 

Since each transformation is orthogonal, the absolute value of the Jacobian in 
each instance is unity. Thus the right member of (2) can now be written as the 
product of s sets of integrals, the sets containing 7,, r2,--- ,7. integrals re- 
spectively. That is, 


P(t, +++ , te) = gilti) +--+ gelts), 


which is equation (1). Hence the theorem. 
Under the conditions of Theorem III, the characteristic function of the 
distribution of 6, is found by direct integration to be 


Tv 


Pv(ty) = (1 — Qio*t,) °. 


6 If the variables in a symmetric quadratic form with matrix A are transformed by a 
linear transformation with matrix B, the new form has the matrix B’ AB. Cf. Bécher, 
p. 129. It should be remarked that these s — 1 successive orthogonal transformations can 
be combined into a single orthogonal transformation with matrix L = L,L2--+ L.1. For 
if, by means of a linear transformation with matrix Z: , we pass from the variables 21, --- , 
ty to the variables z; , --- , y , in which the old variables are expressed explicitly in terms 
of the new, and thence to variables 2; , --- , 2 by means of a linear transformation with 
matrix L, , the transformation with matrix L,L,. will carry us directly from the z’s to the 
z”’s. This extends to any finite number of transformations. Since theproduct of any two 
orthogonal matrices is an orthogonal matrix (and hence the product of a finite number of 
them), we see that the remark is justified. Cf. Bécher, p. 68 and Kowalewski, p.161. Note 
that Bécher expresses the new variables explicitly in terms of the old. 
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f.(6,) de, = dé, = | ete » (t,) dty 


2r Jo 


1 S -3(3) do, 
es € —s 
pe o* o* 
2? r ay 
(3 
so that the variables 0,/o° are distributed in accordance with Chi-square dis- 
tributions with r, degrees of freedom.’ Accordingly, when the conditions of 


Theorem III are satisfied, we may deduce not merely the mutual independence 
of the @, but also the nature of their distributions. 


3. Applications to the analysis of variance. In the analysis of variance, 
N = ab independently observed values of a normally distributed variable are 
classified into a rows and b columns in accordance with some relevant scheme: 

11,712, °** » Lib 
X21 , T22 , oe » L2b 


Lai » Ta2 , oe* »DLab- 


With the notation Z;. , Z., , to denote respectively the arithmetic mean of the 
jth row, the kth column, and the entire set , it is readily seen that 


b 


YE (ee - 2) =b LG. - 2 +aLGs- H+ LD 


1 


(5) 
(xj = Tj. — <5 + z)° 
= A + Oe a 63 


is an identity in the N = ab values of x. It is quite straightforward to exhibit 
each of the three terms in the right member of (5) as a real symmetric quadratic 
form in the N variables z ;, and to show that the ranks are r; = a — 1,7. = b — 1, 


r3 = (a — 1)(b — 1). By the device of adding 6, = + (>> > 2%)? = NZ to 


both members of (5), we have >> >> tik = 0, + 02+ 63+ 6,. Moreover, the rank 
of ists = 1. Thusr, + r2 + 73 + 7% = ab = N and, by Theorem III, we see 
that the four quadratic forms are mutually independent. In particular, 6; , 4 
and @; are independent, and each, meesured in units of o’, is distributed as is 
Chi-square with its appropriate number of degrees of freedom. 


Tue UNIVERSITY oF Iowa. 


7 By the number of degrees of freedom of a real symmetric quadratic form of normally 
and independently distributed variables, we mean the rank of the matrix of the form. 





INDEPENDENCE OF CERTAIN ESTIMATES OF VARIANCE 


REFERENCES 


(1) M. S. Bartierr anp J. WisHart, The distribution of second order moment statistics 
in a normal system. Proceedings of the Cambridge Philosophical Society, 
vol. 28 (1931-32) pp. 455-459. 
The generalized product moment distribution in a normal system. Same journal, 
vol. 29 (1932-33) pp. 260-270. 
(2) W. G. Cocuran, The distribution of quadratic forms in a normal system. Proceedings 
of the Cambridge Philosophical Society, vol. 30 (1933-34) pp. 178-191. 
(3) R. A Fisuer, Applications of Student’s distribution. Metron, vol. 5 (1925) pp. 
90-104. 
Statistical Methods for Research Workers (1934) pp. 210-272. 
(4) S.S. Wixxs, Statistical Inference (1936-37) pp. 38, 44-45. 






























VARIANCE OF A GENERAL MATCHING PROBLEM* 
By JosepH A. GREENWOOD 


Let us match two decks of cards: (A) composed of ¢ distinct groups of s 
identical.symbols each, and (B) a target deck composed of 7; symbols of the 
first kind, 72 of the second, etc., such that 


ytiet--- + = st =n. (1) 
It is not necessary that all the 7z’s be different from zero. 


(a) Forming the Relative Frequency Table. The first part of the paper is 
concerned with forming a 2x2-way table showing the relative frequencies of 


hits and misses of all pairs of cards in the target deck. The notation : indicates 


a miss at the 7th card of the target deck, . ahit. * = j indicates a miss at the 


0 
ith card, with the matching card identical to the jth target card. 
CasEI. ith and jth target cards the same symbol. 


a j Theoretical freq. Weighted freq. 
If 0 then —' 0 n—-s—1l ((¢—1)(n—s—1) 2.1 
0 1 8 (¢—l)s=n—s8s 2.2 
(2) 
1 0 n—s n—s 2.3 
1 1 s—1l s—1 2.4 


Total = t(n — 1) 

z 

But 0 
with a factor (t — 1), giving the last column in (2). 


CasE II. ith and jth target cards different 


occurs in (¢ — 1)/t of the events. Thus we must weight 2.1 and 2.2 








a j. Theoretical freq. Weighted freq. 

If 0 = j then — 0 n—s n—s 3.1 
0=j) 1 s—l s—l1 3.2 
04 j 0 n—s—1 a-a~ 18-8 22 3) 
0Fj 1 s s(t — 2) 3.4 
1 0 n—s—1 n—s—1 3.5 
1 = s 8 3.6 


Total = t(n — 1) 


* Presented to the American Mathematical Society, September 9, 1937. 
1 Read, ‘then out of n — 1 times’. 
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But . = j occurs in 1/(t — 1) of all events =. and : occurs in 1/¢ of all events 


: 4+ . . Therefore entries 3.3 and 3.4 must be weighted with the factor (t — 2), 


and then entries 3.1, 3.2, 3.3 and 3.4 must be weighted with the factor (¢ — 1). 
It is important that the totals of the two parts to be weighted be equal before 
the weighting factors are applied. This gives rise to the last column in table (3). 

Now the number of ways the ith card can be like the jth card of the target 


deck is” 
«a= 2() 


The number of ways they can be unequal is 


i . °,t 
i n 
m= 2 (utd) = (3) - a | (a) 
u<v 
Since the totals of the last columns of the two tables are equal we weight the 
entries of their last columns with a; and a2, respectively. So, combining 3.1, 
3.3 and 3.2, 3.4 we form a; times (2) + ae times (3) to give the new table 


a Z Relative frequencies 
0 0 (n — s — 1)(t — lh + [(¢ — 1)(n — 8 — 1) + Alar 
0 1 (n — s)a, + (n — s — 1)ae 
(5) 
1 0 (n — s)a, + (n — s — l)ae 
1 ] (s — l)ay + Sa . 
Now using the entries from (5) form the 2x2-way table 
Total 
(¢ — 1)(n — s — 1l)ay +} (n — s)ay + (tn —n —t +1) 
[(t — 1)(n — f= 1) (n = _ o 1)ae (ay + ae) 
+1] Qe2 
(n — s)a, + (s — 1)a1 + saz (n — 1)(a1 + az) (6) 
(n — s — l)ae 
(tn —n —t +1) (n — 1)(a, + ae) t(n — 1)(a; + a) 





(a1 + az) 





2 Ifi, < 2 define GS) = 0. 
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(b) Obtaining the Correlation, Variance and Mazimal Conditions. Sub- 


stituting from table (6) into the formulas given by Yule’ for 6 and the coefficient 
of correlation r, we obtain the average correlation 


n 
atone GG) 
(in — mn —t + 1)(on + on) (5) (m —n—0+1) 


n 
2 
tar + (1 — 2) (3) 


(3) em - n—t+1) 
by (4). 


We now give a proof that r is a maximum when 7; = s(j = 1,--- ,#). From 
(7) it is sufficient to show that under the same conditions a is a maximum. 
Let 7; = s + 6;, then 


(7) 


Vent 0. 


a = — (s + 6.)(s + 4) = (;) s+ = bud, by (8). 


u<v 


Assume some 6, ~ 0 and 
1 eee 


7 8 
> 6.6, = 0. 


u<v 


1,-+*,t 1,-*:,¢ 
; s+ Z 6u5y 


u<v 


to both sides of (9). Then 


oe: 2 1,--°,¢ 1-**,t 
(= i.) >Le+ DL ws 


u<v 


or 0 2 a positive number. This necessarily implies the desired result. 


*Yule, G. U. An Introduction to the Theory of Statistics, London: Griffin and Co., 
1927, pp. 216-217. The table can be symbolized with 
Total 


6 be — (c2b3/cs) 


He then gives r = 8c3/4/c,C2asbs, 


the correlation coefficient. 
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Yule‘ gives an expression for the variance in a situation which includes the 
present problem as a special case, to be 


o = npg[l + r(n — 1)], (11) 


where r is the average correlation between all pairs of variables. Substituting 
our result (7) in (11) with p = 1/t gives the desired variance. 

It is interesting to note that when 7; = s, (j = 1, --- , 4) rreduces to 1/(n — 1)* 
giving 


2_ m(t—1) _ n o 
#(n—1) n—-1’ 


2. ‘ ‘ 
where o; is the variance of the binomial case.° 


DvuKE UNIVERSITY. 


‘Op. cit., p. 286. 

5 Concerning this special case see also Bartlett, M. 8. Properties of sufficiency and 
statistical tests. Proc. Royal Soc. A. 1937, CLX, 268-282. 

Olds, E.G. A moment-generating function useful in certain matching problems. Abstract 
No. 428, Bull. Amer. Math. Soc. 1937, XLIII, 779. 





THE LARGE-SAMPLE DISTRIBUTION OF THE LIKELIHOOD RATIO 
FOR TESTING COMPOSITE HYPOTHESES’ 


By S. 8S. Wi1LkKs 


By applying the principle of maximum likelihood, J. Neyman and E. §. 
Pearson’ have suggested a method for obtaining functions of observations for 
testing what are called composite statistical hypotheses, or simply composite 
hypotheses. The procedure is essentially as follows: A population K is assumed 
in which a variate x (x may be a vector with each component representing a 
variate) has a distribution function f(z, 6; , 62, --- 0,), which depends on the 
parameters @,, 02--- 6,. A simple hypothesis is one in which the 6’s have 
specified values. A set © of admissible hypotheses is considered which consists 
of a set of simple hypotheses. Geometrically, 2 may be represented as a 
region in the h-dimensional space of the 6’s. <A set w of simple hypotheses is 
’ specified by taking all simple hypotheses of the set 2 for which 6; = @;, 7 = 
m+i1,m+2,---h. 

A random sample O, of n individuals is considered from K. O, may be 
geometrically represented as a point in an n-dimensional space of the x’s. The 
probability density function associated with O, is 


(1) P = TT flee, 015015 +++ 6) 


Let Po(O,) be the least upper bound of P for the simple hypotheses in Q, and 
P..(O,) the least upper bound of P for those inw. Then 


_ P.(On) 
(2) : = B500,) 


is defined as the likelihood ratio for testing the composite hypothesis H that 
O, is from a population with a distribution characterized by values of the 6; 
for some simple hypothesis in the set w. When we say that H is true, we shall 
mean that O, is from some population of the set just described. In most of the 
cases of any practical importance, P and its first and second derivatives with 
respect to the 6; are continuous functions of the @; almost everywhere in a certain 
region of the 6-space for almost all possible samples O,. We shall only consider 
the case in which P»(O,) and P.(O,) can be determined from the first and 
second order derivatives with respect to the @’s. 


1 Presented to the American Mathmatical Society, March 26, 1937. 
?Phil. Trans. Roy. Soc. London, Ser. A, Vol. 231, p. 295. 
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A considerable number of currently used statistical functions for making tests 
of significance can be expressed in terms of \ ratios, and in many cases involving 
normal distribution theory, the exact sampling distribution of \ is known. 
However, it is often useful when dealing with large samples to have an approxi- 
mation to the distribution of A. We shall consider such an approximation for 
those cases (which include most of the ones of any practical importance) in 
which optimum estimates of the 6’s exist. That is, we shall assume the existence 
of functions 6,(2; , --- z,)(maximum likelihood estimates of the @;) such that’ 
their distribution is 


| cj (* a] 5 C4 j2iej 


(3) (2n)*2 = (1 + ¢) dz, eee dzp, 


2 
where 2; = (6; — 0:)-/n, cj = —E (5 was), E denoting mathematical expecta- 
4 OU; 


tion, and ¢ is of order 1/+/n and ||c;;|| is positive definite. Denoting (3) by 
Jdz,dz, --- dz, , and differentiating J with respect to 6 , we get 


Uf 4. ted. Stews vB) b= 
(4) Ate 30, by am aaj + Vn De cuit; J, k=1,2, h 
Since c;; = O(1) and |c;;| # 0, it can be seen from (4) that the values of 6 
which maximize J differ from 6,, k = 1, 2,---h, by terms of order 1/V/n. 
Therefore, the maximum P2(O,) of J with respect to the 4 is oon (1 + 4’, 
where ¢’ = O(1/+/n). 

To get P.(O,), we let 0; = 0;,71 = m-+1,m 4+ 2, --- h, and note that J can 
be written as 


| Coij ! 7 ; B csatieang ’ 
(5) Jo= (nya " (1 + ¢o) 


where 
h 


(6) xo = 2 Ci52i2i $0 = O(1//n) 
+,7=m 

and || c;; || is the inverse of the matrix obtained by deleting the first m rows and 

first m columns from || ¢;; |{"' and z; = zi — L,;, L; being a linear function of 

90,m+1 °*+ On, and ¢o;; is the value of c,; with 6; = 0;,7 =m+1,m+2,.---h, 

that is, when H is true. Taking the maximum P,(O,) of expression (5) with 

respect to 6; , 62, --+ Om, we get 


y(* 2 ” ” 
(7 Py = (atthe eMC + 40) $0 = O(1/-Vn) 


*For conditions under which the 6’s exist which are distributed according to (3), see 
J. L. Doob, Probability and Statistics, Trans. Amer. Math. Soc. Vol. 36, p. 759-775. 
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Hence, when H is true, we have, from (5) and (7) 


(8) ne coe = (1 + O(1/+/n)). 


Therefore, except for terms of order 1/+/n, 
(9) — 2logd = xo. 


Now, the characteristic function of —2 log } is 


¢(t) aes E(e"* log ))) oe / = [ aoe siroerve ae, 7 dz, 
10 
, : | cs |* = es ja52;+x9 (it) 
= Gry / ne / ew (1 + O(1/Vn)) des «++ dep. 


It can be shown that on any finite interval | t | < a, $(¢) approaches uniformly, 
as n— o, the function 


(11) QFa-w*. 


But (11) is the characteristic function of any quantity distributed like x’ with 
h — m degrees of freedom. 

We can summarize in the 

Theorem: If a population with a variate x is distributed according to the probabil- 
ity function f(x, 0; , 02 --+ On), such that optimum estimates 6; of the 0; exist which 
are distributed in large samples according to (3), then when the hypothesis H is 
true that 0; = %:,71 = m-+1,m-+ 2, --- h, the distribution of — 2 log X, where d 
is given by (2) is, except for terms of order 1/+/n, distributed like x* with h — m 
degrees of freedom. 


Princ&ton UNIVERSITY, 
PRINCETON, N. J. 





ON DIFFERENTIAL OPERATORS DEVELOPED BY O’TOOLE 
By M. Z1aup-DIN 


1. O’Toole in his paper ‘Symmetric Functions and Symmetric Functions of 
Symmetric Functions’ [Ann. Statist. 2. (1931)102-49], has expressed Monomial 
Symmetric Functions }>?1???* ... , in terms of power-sums, s, . 

The Monomial Symmetric Functions can be written in partition notation as 
#1 ia *3...) where ki, ke, --- denote the repetitions of parts. 

To express (‘12,3 ...) as a function of s,, O’Toole has developed operators 
d, and D, , connected by the formulae,’ 


d 

ds,’ 

rd, = (HD (HD & = lr Da De --- 
kilkal => 


rip, = 2tidadi «+ 
ia: 


where kA +. &B+--- =r 
ky + ke +.--=k. 


In this paper it will be shown that these operational relations are easily 
deduced from the operators d, and D, of Hammond, used for expressing Mono- 
mial Symmetric Functions as functions of Elementary Symmetric Functions, a, . 

For the sake of distinction I shall use g, and Q, for the operators employed by 
O’Toole and keep d, and D, for Hammond’s Operators. 

Macmahon has dealt with Hammond’s operators in his Combinatory Analysis 
Vol. I. Cambridge University Press (1915), where they are defined’ as 


d, = 


_ _d d d 
D, = ~ (di) and o? +e +69? -* (1). 


2. It is known’ that 


log (1 — aiz + mx” —a;2°+.---) = —(s2+ 5 os" + tee +i a2' + .), 


1O’Toole, Loc. cit., p. 120. 
2? Macmahon. Comb. Analysis. I. 27-28. 
* Tbid., p. 6. 








64 M. ZIAUD-DIN 













Now operate on the right hand side with d, , and with its equivalent th (1) on 
the left hand side. Equating coefficients of zx’ on both sides, we obtain, 
° d,8 = (—1)"'r; d,s, = 0 when r # k, 

which yields 


ee ee x tunis. (2), 








The operator Q, exactly behaves like D,. From the formula‘ 
d, faa D,d,-1 + Ded,-2 — Ds3d,—~s + ee + (—1)'rd, = 0, 


which is in complete correspondence with Newton’s recurrence relation, we 
derive 









dy = 
ds 


D, 
pi — 2D (C) 
Di-— 3D, De + BDo. 


Poe eee UES ESE COSECOCSOCOOOSOSOOCOOOOOOOCOS OE 






By multinomial theorem 
g = SMH = VD De --- 
. kylke! -- 
using (2) we at once get 
rg, = (— Itt SLED = Qe GF 
: ky!ke! --- 
which is the result (A) obtained by O’Toole. 
From (C), D, follows in terms of d, and thence with (2) Q, can be expressed in 
terms of g,. Using multinomial theorem we arrive at 


Iq’ Ke 
r! . 
r!Q, Zz A B 


é ky! ke! --- 
which is (B). Hence both the results of O’Toole have been deduced. 


















3. In his second paper’ O’Toole defines symmetric functions for more than 
one system of Variates. I call such symmetric functions Hyper Symmetric 
Functions. 

The Hyper operators are developed to express Hyper symmetric functions in 
terms of hyper power-sums. They are defined by O’Toole by the following 
relations, taking into consideration two systems of variates only, 

d’ 


d,, = —; 
Pq ’ 
ASpq 







4p. 29. 
6 Ann. Stat. 3. (1932), 56-63. 
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d., = i, (— 1k Pz 1)!Diia, i... Sides 
= klk! --- 
i. roa 


Doo ee 


ky! ke! - 
where kip; + kepe + --- p 
kigi + kee + --- = Q 


These relations readily follow from Macmahon’s’ hyper operators gp, and 
Gpq. These operators came into existence with the problem of expressing 
hyper symmetric functions in terms of hyper elementary symmetric functions 
and they are connected by the following relations. 


L (epee @t9-=D! _ V(-1)%k—D! on oe 
. "plat 9 ~ ~~ i.. ~orm * 


(- "4, 


To Limt+a- 1 (mate —I (-p 


k k 
pilg:! pilgal | Rilke === Ime) Orne) * +-- 
Macmahon’ has also shown that 


yes pq! . 
@+¢- 1D! 


JpqgSpq = (-1 


from which we get 


ae ce p'q! 
Gog = ( 1) (p 4 q- 11 (3) 


The operator G behaves like D of O’Toole. . Now using (3) we derive from (I) 
the result (A’) arrived at by O’Toole without reference to Macmahon. Simi- 
larly from II. using (3) (B’) is deduced. 


UNIVERSITY COLLEGE, SWANSEA, WALES. 


* Macmahon. Comb. Analysis. Vol. II. Cambridge University Press (1916), p. 302. 
7 Macmahon Op. Cit., p. 304. 








































GRADUATION BY A TRUNCATED NORMAL 


By NatTHAN KEyYFITz 


Below is a table for finding the constants of a truncated normal by the equa- 
tion of moments. Karl Pearson* gives such a table for the case in which the 
data are to be fitted to the “tail” (i.e. less than half) of a normal curve but I do 
not believe that the formulae for a distribution consisting of more than half of a 
normal curve have before been tabulated. 

The table below was calculated primarily for an investigation being carried 
out on the duration of unemployment. The Canadian Census of 1931 reported 
the number of persons losing 1-4, 5-8, - - - 49-52 weeks in the course of the year 
June Ist, 1930 to June Ist, 1931, by various classifications, (industry, province, 
age, etc.). 

The tendency to report even numbers of months on the part of the enumerated 
population was evident in the result, and some kind of graduation was necessary 
for an interpretation. After some experiment a part of a normal curve was 
settled upon as the simplest and generally most satisfactory representation. 

It was found that among the classes of workers in which unemployment is 
high the curve is more advanced ,—i.e. the mode is at a higher number of weeks,— 
than in the classes where unemployment is low. In many cases, (in most 
groupings of female workers for instance) where unemployment is relatively very 
low the modal point of the uncurtailed normal stands at a negative number of 
weeks,—for these cases the fitting is to a true tail and the tables of the Biometric 
Laboratory were used. 

Details of the results of the investigation will be published shortly in the 
Unemployment Monograph of the Dominion Bureau of Statistics. Meanwhile, 
this table will be of use as the complement of Pearson’s tabulation which is only 
suitable for ¥, = .5708. 


Table for finding the constants of a truncated normal by the equation of moments 


x vi _ oh al. __ Ave 

0 .5708 — .0180 1.2533 — .0562 
a 5528 ‘— 0183 1.1971 — .0543 
3 5345 — .0188 1.1428 — .0526 
3 .5157 — .0190 1.0902 — .0506 
4 .4967 — .0193 1.0396 — .0487 
5 .4774 — .0195 .9909 — .0467 





* Tables for Statisticians and Biometricians, page 25. 
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A v1 _ ahr ve einilliais 
6 .4579 — .0195 .9442 — .0449 
= .4384 — .0196 .8993 — .0428 
8 .4188 — .0196 .8565 — .0409 
9 .3992 — .0194 .8156 — .0390 
1.0 .3198 — .0192 .7766 — .0370 
Bad . 3606 — .0189 . 7396 — .0351 
1.2 3417 — .0185 . 7045 — .0332 
1.3 3232 — .0180 .6713 — .0315 
1.4 .3052 — .0175 .6398 — .0296 
1.5 2877 — .0170 .6102 — .0279 
1.6 . 2707 — .0163 . 5823 — .0263 
1.7 . 2544 — .0156 . 5560 — .0246 
1.8 . 2388 — .0148 .5314 — .0232 
1.9 . 2240 — .0141 . 5082 — .0216 
2.0 . 2099 — .0134 .4866 — .0204 
2.1 . 1965 — .0126 .4662 — .0190 
2.2 . 1839 — .0118 .4472 — .0178 
2.3 .1721 — .0110 .4294 — .0166 
2.4 .1611 — .0103 .4128 — .0156 
2.5 . 1508 — .0096 .3972 — .0146 
2.6 .1412 — .0089 .3826 — .0137 
2.4 . 1323 — .0083 . 3689 — .0128 
2.8 . 1240 — .0076 .3561 — .0120 
2.9 . 1164 — .0071 3441 — .0113 
3.0 . 1093 .3328 
3.5 .0813 . 2856 
4.0 .06246 . 24999 
4.5 .049379 . 222221 
5.0 .0399997 . . 1999999 


Let d = distance of centroid of actual distribution from: point of truncation. 
Let = = standard deviation of distribution about its mean. Then y, = a , 


Hence corresponding x’ and ¥2 may be found. 

Then o = dye, where o = standard deviation of uncurtailed normal. 

And x = x’o, where x = origin of uncurtailed normal. 

N.B. The point of truncation is taken for the origin in the original distribution. 


SocraL ANALYSIS BRANCH, DOMINION BUREAU OF STATISTICS, 
Orrawa, CANADA 











REPORT OF THE ANNUAL MEETING OF THE INSTITUTE OF 
MATHEMATICAL STATISTICS 


The annual meeting of the Institute of Mathematical Statistics was held on 
Wednesday and Thursday, December 29-30, 1937, in Indianapolis, Indiana, in 7 
conjunction with the meetings of the American Mathematical Society and | 
associated organizations. 

The Wednesday morning session was devoted to applications of statistics to 
industry and engineering. On Thursday morning, the Institute held a joint 
session with the Mathematical Society for the presentation of voluntary papers 
on probability and statistics. This session was immediately followed by 
another of the Institute for two invited addresses. These addresses were 
“The theory of general means” by Professor E. L. Dodd, and “On the inde- 
pendence of certain estimates of variance” by Professor A. T. Craig. Professor 
P. R. Rider was in charge of arranging the program. 

On Thursday noon, there was a luncheon at the Marott Hotel for members of 
the Institute and their guests. After the luncheon, Professor H. L. Rietz spoke 
on “The future of the Institute in relation to mathematical statistics.” 

At the business meeting, which followed the Wednesday morning session, 
President Shewhart announced that these officers had been elected for 1938: 
President, B. H. Camp, Wesleyan University; Vice-Presidents, P. R. Rider, 
Washington University, and 8. S. Wilks, Princeton University; Secretary- 
Treasurer, A. T. Craig, University of Iowa. The Institute voted to hold its 
1938 meeting with the American Statistical Association. The meeting will be 
in Detroit, Michigan, in December of this year. 

ALLEN T. Crate, Secretary. 





